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Lidar  Measurement  of  Optical  Turbulence: 
Theory  of  the  Crossed  Path  Technique 


1.  INTRODUCTION 

With  the  advent  of  lasers  and  their  widespread  application,  the  importance  of  atmospheric 
optical  turbulence  has  grown.  Atmospheric  turbulence  degrades  the  propagation  of  laser  beams 
through  the  atmosphere.  Specifically,  wind  fluctuations  induce  density  or  temperature  fluctuations. 
These  temperature  fluctuations  in  turn  result  in  fluctuations  of  the  atmospheric  refractive  index  on 
scale  sizes  from  a  centimeter  (or  less)  to  tens  of  meters.  Thus,  the  refractive  index  is  not 
homogeneous  across  a  laser  beam  wavefront  larger  than  approximately  a  centimeter.  The  resulting 
effects  of  the  turbulence  are  a  loss  of  coherence  in  the  beam,  beam  broadening,  beam  wander,  in¬ 
tensity  fluctuations  or  scintillations  and  intensity  hot  spots  within  the  beam. 

Analogous  effects  occur  in  passive  imaging  systems.  Mitigation  of  these  effects  in  lasers  and 
imaging  is  the  goal  of  adaptive  optics  compensation  methods.  Such  methods  utilize  deformable 
mirrors  to  compensate  for  the  turbulence  induced  phase  distortions.  The  performance  and  charac¬ 
terization  of  adaptive  optics  systems  requires  a  knowledge  of  the  fundamental  turbulence  parame¬ 
ter,  C„^,  as  well  as  auxiliary  meteorological  and  wind  data.  These  parameters  allow  specification 
of  the  number  of  mirror  actuators,  the  maximum  allowable  look-ahead  or  lag  angle,  and  the  band¬ 
width  requirements  (or  maximum  lag  time)  of  the  adaptive  system.  A  point  that  requires  emphasis 
is  that  these  systems  require  the  characterization  of  an  optical  path.  Such  paths  are  typically  of  the 
order  of  a  meter  in  diameter.  Furthermore,  such  systems  require  the  characterization  of  the  optical 
path  over  short  time  scales.  Of  course,  these  paths  need  to  be  specified  in  terms  of  ensemble  aver¬ 
ages;  the  specification  of  an  instantaneous  wavefront  fluctuation  is  required  only  for  the  entire 
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propagation  path  and  provides  the  specific  information  for  performing  phase  conjugation  using 
adaptive  compo  lents. 

The  requirements  for  the  specification  of  the  atmospheric  effects  are  in  terms  of  small 
volume  averages.  Such  high  resolution  measurement  ctqpability  is  also  needed  for  more  fundamen¬ 
tal  studies  of  atmospheric  turbulence,  specifically,  the  fine  scale  spatial  and  temporal  structure  of 
turbulence.  These  issues  of  turbulence  structure  pertain  to  stationarity,  intermittency  and  the  as¬ 
sumption  of  frozen  turbulence.  One  measurement  system  that  seems  (at  Hrst  consideration)  to 
have  the  potential  for  providing  such  high  resolution  capability  involves  the  use  of  lasers  for  atmo¬ 
spheric  remote  sensing,  that  is,  lidar.  Since  the  early  days  of  lidar,  atten^ts  have  been  made  to  use 
lidar  to  measure  the  atmospheric  refractive  index  fluctuations.  Very  limited  success  has  been 
achieved. 

In  this  report,  the  various  existing  methods  for  measuring  optical  turbulence  and  its  parame¬ 
ters  are  surveyed.  This  demonstrates  the  source  of  our  existing  knowledge  of  atmospheric  optical 
turbulence  and  outlines  the  limitations  of  the  systems  and  various  measurement  issues.  A  second 
brief  survey  of  the  various  approaches  to  utilizing  lidar  for  profiling  follows.  The  last  several 
sections  are  devoted  to  the  proposal  of  a  concept  for  lidar  turbulence  sensing. 


2.  METHODS  OF  MEASUREMENT  OF  OPTICAL  TURBULENCE 

Measurements  of  optical  turbulence  have  been  made  for  several  years  via  radars,  in  situ  in¬ 
strumentation  carried  by  balloons,  and  passive  optical  instrumentation.  All  the  existing  methods 
are  limited  in  various  ways. 

Radars  suffer  from  high  cost,  and  the  large  space  requirements  of  the  antennas  eliminate 
portability.  The  maximum  vertical  resolution  is  of  the  order  of  150  m  in  the  latest  generation  Dop¬ 
pler  systems',  but  horizontal  resolution  is  dictated  by  beam  divergence  and  consequently  horizon¬ 
tal  resolution  is  altitude  dependent.  Divergence  is  usually  given  in  degrees  or  tens  of  degrees;  the 
Flatland  VHF  radar^  has  a  divergence  of  3.6  degrees  or  626  m  at  10  km  while  the  Sousy  VHF 
radar^  has  8  degrees,  or  1.4  km  horizontal  resolution  at  10  km.  Additionally,  radars  typically 
measure  a  weak  signal  and  require  long  averaging  times;  data  is  usually  reported  as  15  minute  or  1 
hour  averages.  Shorter  averaging  times  are  possible  by  using  larger  vertical  range  cells  that  are  as 
large  as  1  km  in  some  cases.  Radars  measure  the  average  turbulence  in  a  large  volume,  many 
times  larger  than  a  beam  or  imaging  line  of  sight.  Radars  also  suffer  from  two  other  critical  flaws. 
First,  they  cannot  measure  within  the  first  few  kilometers  above  the  ground.  This  region  includes 
the  atmospheric  boundary  layer  which  has  a  dominant  contribution  to  various  optical  effects  such 
as  beam  wander,  broadening  and  loss  of  coherence.  Second,  radars  are  typically  limited  to  a  maxi- 


'Hocking,  W.K.  (1986)  Ob*erv»iion  «nd  measurement  of  turbulence  in  the  middle  atmosphere  with  a  VHF  radar,  J.  Amos,  and  Terns. 
Physics,  411:655-670. 


^reen,  J.L.,  Beland,  R.R.,  Brown,  J.H.,  Clark,  W.L.,  Eaton,  F.D.,  Favier,  L.D.,  Gage,  K.S.,  Hatch,  W.H.,  Hines,  J.R.,  Murphy.  E.A..  Nastrom, 
O.D.,  Peterson.  W.A.,  VanZandt,  T.E.,  and  Wamock,  J.M.  (1989)  Comparisons  of  Refractivity  Measurements  from  the  Flatlands  VHF  Radar 
with  Other  Measurement  Techniques,  IVoceeedtngs  of  the  24*'  Conference  on  Radar  Meteorology,  Tatlaha.ssee,  FL,  March  27-31. 

^Eaton,  F.D  ,  Peterson,  W.A.,  Hines,  J.R.,  Peterman,  K.R..  Good,  R.E.,  Beland,  R.R.  and  Brown,  J.H.  ( 1988)  Compariaona  of  VHF  radar,  opti¬ 
cal  and  temperature  fluctuation  measurements  of  r^,,  and  9^,  Thtor.  Appl.  Climatol.,  39: 17-29. 
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mum  range  of  20  km,  but  this  is  a  less  serious  flaw  since  turbulence  is  of  less  importance  in  the 
stratosphere  in  most  applications.  This  range  limitation  can  be  increased  by  increasing  the  radar 
power.  Finally,  the  physical  principle  in  turbulence  measurement  by  radars  is  the  backscattering  of 
the  beam  by  inhomogeneities  in  the  radar  refractive  index  that  are  the  same  scale  as  the 
wavelength  (a  VHF  radar  operates  typically  at  50  MHz,  or  A,  =  6  m,  which  is  assumed  to  be  within 
the  inertial  range  of  turbulence).  Note  that  the  backscattering  results  from  refractive  index  fluctua¬ 
tions  at  radar  wavelengths,  not  optical/IR  wavelengths.  For  radar  wavelengths,  the  refractive  index 
fluctuations  are  dominated  by  humidity  fluctuations,  while  for  optical/IR  wavelengths,  the 
temperature  fluctuations  are  dominant.  To  convert  from  radar  to  optical  turbulence  requires  the 
measurements  of  pressure,  temperature  and  humidity.  The  traditional  approach  has  been  to  use 
nearby  rawinsonde  profiles.  These  are  low  resolution  (about  1  km  vertically)  and  are  not  in  the 
same  volume.  This  is  a  serious  drawback  if  the  goal  is  to  use  radar  to  characterize  an  optical  path. 

Balloon-borne  measurements  of  optical  turbulence  have  been  made  for  some  time^.  The 
prevalent  technique  involves  the  measurement  of  temperature  fluctuations.  Since  optical  refractive 
index  fluctuations  arise  entirely  from  temperature  fluctuations,  these  methods  provide  a  direct 
measurement  of  Cn^.  A  further  advantage  of  these  measurements  is  their  ability  to  simultaneously 
measure  meteorological  and  wind  data  through  the  inclusion  of  a  standard  radiosonde  in  the  in- 
stmment  package.  The  vertical  resolution  of  these  methods  is  determined  by  statistical  estimation 
considerations  and  is  typically  of  the  order  of  50  meters.  A  major  problem  with  these  in-situ  meth¬ 
ods  is  that  they  are  not  made  along  any  optical  path,  but  rather  a  path  determined  by  the  wind.  A 
second  drawback  is  that,  although  they  provide  the  highest  spatial  vertical  resolution,  balloon- 
borne  measurements  also  provide  the  poorest  temporal  resolution;  a  balloon  flight  takes  on  the  or¬ 
der  of  an  hour.  Finally,  daytime  measurements  in  the  stratosphere  and  upper  troposphere  present 
problems-'  that  are  likely  attributable  to  solar  heating  of  the  temperature  probes. 

Most  passive  optical  measurement  techniques  utilize  stars  as  sources  and  measure  some  tur¬ 
bulence  effect  on  the  propagation  of  the  starlight.  Thu.s,  a  stellar  scintillometer  measures  the  star’s 
intensity  fluctuations  in  a  telescope  aperture  to  infer  an  isoplanometer  measures  isoplanatic 
angle  from  a  star;  and  an  MTF  device  measures  the  transverse  coherence  length,  r^,  by  measuring 
the  blur  circle  of  a  star’s  image  in  a  telescope.  These  are  the  most  widespread  techniques  for 
measuring  the  vertical  stmeture  of  turbulence.  In  principle,  they  are  capable  of  performing  tur¬ 
bulence  measurements  along  an  optical  path  provided  there  is  a  conveniently  located  star.  All 
these  techniques  are  capable  of  high  temporal  resolution,  limited  only  by  the  statistical  estimation 
requirements.  However,  both  the  MTF  and  isoplanometer  measure  a  moment  or  integral  of  C„^; 
they  provide  no  information  about  the  distribution  or  layering  of  turbulence.  Thus,  their  altitude 
resolution  is  nonexistent.  The  scintillometer  does  provide  some  altitude-resolved  information,  but 
the  altitude  resolution  is  poor,  of  the  order  of  a  few  kilometers.  Further,  scintUlometers  tend  to  re¬ 
quire  averaging  times  of  the  order  of  5  or  more  minutes.  A  serious  drawback  of  the  scintillometer 
is  its  inability  to  measure  turbulence  in  the  first  kilometer,  a  very  important  region  as  described 
earlier.  Another  limitation  of  the  scintillometer  is  its  restriction  to  night  operation.  Both  the 
isoplanometer  and  the  MTF  devices  have  been  adapted  for  daytime  use  as  well. 


^rown,  J.H.,  Good.  R.E.,  Bench,  P.M.  and  Faucher,  O.  (1982)  Sondr  Measiiremenis  for  Comparative  Measurements  of  Optical  Turbulence. 
Air  Force  Geophysics  Lab.,  APOL-TR-82'0079,  ADA  1 18740.  NTIS. 


■'WaHers,  D.L.,  Weilekamp,  M.R.,  Beland.  R.R.,  Good,  R.E.  arid  Murphy,  E.A.  (1990)  Single  fVobe  Optical  Turbulence  Profile  Measurement 
System,  Symposium  on  Propagation  Through  Turbulence,  Annual  Meeting  of  the  Opt.  Soc.  Am.,  Bostoit,  MA,  November  4-9. 
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This  survey  of  various  measurement  techniques  is  intended  to  present  the  capabilities  and 
limitations  of  the  various  approaches.  It  emphasizes  the  methods  currently  in  routine  use  for  atmo¬ 
spheric  characterization.  The  survey  is  not  complete,  as  there  are  other  methods  that  have  been 
used.  One  proposed  optical  technique  uses  the  scintillation  from  a  satellite-based  light  source  for 
C„^  profiling*.  Another  optical  technique  that  has  been  presented  in  the  literature  utilizes  scintilla¬ 
tion  from  binary  stars  to  obtain  C„^  profiles  with  up  to  1  km  vertical  resolution’.  Due  to  the  rarity 
of  binary  stars,  this  technique  is  limited  and  cannot  be  used  to  measure  along  a  specific  optical 
path.  A  variation  of  this  technique  will  be  discussed  in  more  detail  subsequently  in  the  context  of 
lidar  methods.  The  purpose  of  this  survey  has  been  to  demonstrate  the  capabilities  and  limitations 
of  existing  methods  and  to  show  that  there  is  currently  no  method  capable  of  measuring  along 
a  specific  optical  path  with  sufficient  spatial  and  temporal  resolution. 


3.  LIDAR  APPROACH  TO  MEASUREMENT  OF 


Given  the  requirement  of  characterizing  a  specific  optical  path,  the  use  of  lidar  is  an  obvious 
candidate.  Lidar  methods  fall  under  the  category  of  active  optical  techniques.  The  generic  lidar  ap¬ 
proach  to  any  atmospheric  measurement  is  to  use  a  laser  pulse  that  is  range  gated  at  the  receiver  so 
that  the  time  of  arrival  of  the  pulse  at  the  receiver  determines  the  altitude  of  the  backscatterers.  A 
one  microsecond  pulse  gating  is  typical  of  lidars  and  translates  into  a  150  m  vertical  resolution.  In 
effect,  this  receiver  gating  allows  us  to  treat  the  lidar  as  an  array  of  sources  separated  by  150  m 
along  the  optical  path. 

A  general  formulation  of  a  lidar  approach  to  C„^  measurement  will  be  presented.  A  lidar  sys¬ 
tem  will  measure  some  optical  effect,  say  g,  from  the  gated  return  at  altitude  z,.  In  general,  the  ef¬ 
fect  g  will  be  determined  by  the  entire  optical  path  from  the  transmitter  (z  =  0)  to  the  backscatterer 
(z  =  Zj)  and  then  back  to  the  receiver  (z  =  0).  TTie  optical  effect  can  be  written  as  an  integral  along 
the  optical  path: 


g(Zi)  = 


r'c„’(h)W(X,Zi,h)dh, 

0 


(1) 


where  W(A,,Zj,h)  is  the  particular  optical  effect’s  path  weighting  function  at  wavelength  X,,  ob¬ 
served  at  z  =  0  and  backscattered  at  Zj.  The  weighting  function  measures  the  contribution  to  this 
effect  arising  from  optical  turbulence  at  z  =  h.  We  can  write  this  down  as  a  sum  of  integrals  over 
the  range  gates  as: 

i  z- 

g(Zi)  =  Z  'c„’(h)W(X„Zi,h)dh 
j=l  Zj-, 


.Zj.h)  dh, 


(2) 


*Kr»use-Potstoiff,  J.  bikJ  Wulters,  D.L.  (1990)  Refractive  turbulence  profiling  uaing  an  orbiting  light  aource,  Appl.  Opt.,  29:1877-1885. 

^Azouit,  M.,  Vemin,  J.,  Barletti,  R.,  Ceppatelli,  O.,  Righini,  A.  and  Speroni.  N.  (1980)  Remote  senaing  of  atmoapheric  turbulence  by  meana  of 
a  faal  optical  method;  A  compariaon  with  aimullaneoua  in  situ  meaaurementa,  J.  Appl.  Meteor.,  19:8.54-838. 
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where  C„^(Zj)  is  the  weighted-average  C„^  over  the  altitude  bin  (Zj.|^j)  defined  by 


J^'c„2(h)W(>.,Zi4i)dh 

-ft] _ 


(3) 


dh 


Strictly,  we  should  write  C„^(Zj,Z})  to  show  the  dependence  on  Zj.  However,  the  dependence  of  the 
weighted  average  on  Z;  is  assumed  to  be  weak  over  the  range  (Zj.,,Zj).  Additionally,  if  C„^  is  con¬ 
stant  over  (Zj.,,Zj),  the  dependence  disappears  and  it  will  be  assumed  that  C„^  is  constant  over 
(Zj.,,Zj).  This  last  form  allows  us  to  write  the  result  in  the  suggestive  matrix  form: 
i 


g(Zi)  =  I  Wi^CCzp, 
j=l 


(4) 


or  in  terms  of  vectors  G  and  C  and  matrix  W: 


G  =  WC, 


(5) 


where  the  matrix  W  has  elements  Wjj  defined  by 


W(X,Zi,h)  dh. 


(6) 


and  the  vector  elements  of  G  and  of  C  are  g(Zi)  and  C„^(Zj),  respectively.  The  solution  giving 
C„^(Zj)  in  terms  of  the  measured  quantities  g(Zj)  is  then  formally  given  by: 


C  =  W  G  . 


(7) 


The  formulation  is  readily  extended  to  included  measurement  errors  in  terms  of  an  error  vector  e. 
This  yields  a  formal  least  squares  type  problem  whose  solution  properties  are  the  critical  question. 
The  numerical  stability  and  problem  of  ill-conditioning  are  of  central  importance.  In  general,  the 
matrix  W  has  elements  that  satisfy 


Wjj  =  0  for  j  >  i. 


(8) 


This  simply  states  that  turbulence  at  altitudes  higher  than  the  backscatter  volume  has  no  effect.  All 
the  upper  half  matrix  elements  are  zero  and  the  formulation  gives  a  triangular  system  of  equations. 
The  formal  solution  of  this  system  of  equations  is  given  by* 


i-1 


C(Zi>  =  (W„)-’[g(Zi)  - 1  W,J  O))]-  (9) 

j=l 

There  exist  two  broad  approaches  to  turbulence  measurements  via  lidar.  The  two  cases  are 
where  the  optical  effect  arises  from  phase  distortions  and  where  the  effects  arises  from  amplitude 


'Oahlquist,  O.  and  Bjorck,  A.  {\914) NHinerical Methods,  Prenlice-Hall,  Englewood  Cliffs,  NJ. 
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distortions;  these  are  considered  in  turn.  To  simplify  the  discussion,  it  is  assumed  that  the  back- 
scattering  volume  has  a  cross-sectional  area  (tranverse  to  the  propagation  direction)  that  is  suffi¬ 
ciently  small  so  that  it  can  be  considered  a  point  source.  The  complexity  of  finite  backscatter  spot 
size,  although  important  in  describing  the  optical  effect,  is  not  impoitant  in  this  treatment,  as  we 
are  presenting  "first  order"  considerations.  The  second  major  simplifying  assumption  is  that  there 
are  no  turbulence  effects  on  the  uplink;  that  is,  turbulence  only  effects  the  downlink  propagation 
from  the  backscatterers  to  the  receiver.  This  is  clearly  uruealistic,  but  again,  we  wish  to  discuss  is¬ 
sues  associated  with  a  "first  order"  reproach.  The  effect  of  these  two  assumptions  is  that  we  are 
addressing  the  following  problem:  we  have  a  linear  arrray  of  point  sources  separated  vertically  by 
a  distance  A,  and  we  seek  to  determine  C„^  vertically  from  this  array. 

Consider  first  a  phase  process,  chosen  to  be  image  dancing  for  concreteness.  For  a  point 
source  at  z  observed  by  a  receiver  located  at  0  and  of  diameter  D  ,  the  mean  square  wander,  <5^>, 
is  given  by 


<52(z)>  =  2.91  f^D  *'^  f  C„2(t|)  [i-(ti/z)]5'3  dti. 

*^0 


(10) 


The  expression  gives  the  wander  measured  in  the  focal  plane  of  the  receiver  of  focal  length  f,  pro¬ 
vided  the  receiver  diameter  is  within  the  inertial  range.  In  our  formulation,  we  have  that 


g(Zi)  =  <62(Zi)>,  (11) 

Wij  =  J  '  W(X,Zi,h)dTi 

=  2.91  f  D  f  [i-(ti/z)]5/-^  dt;.  (12) 

This  can  be  evaluated  explicitly  to  give 

Wj.  =  2.91  f  D-'/^  (3z/8)  {[l-(Zj.,/Zi)]»/^ .  (13) 

A  case  of  sp)ecial  interest  is  where  J«i,  that  is,  the  altitude  range  bin  is  far  from  the  backscatter 
region.  An  expansion  yields 

Wij »  2.91  f  (Zj  -  Zj.,)  =  2.91  f^  A,  (14) 


where  A  =  (zj  -  Zj.,)  is  the  width  of  the  range  gate.  Note  that  the  matrix  elements  far  from  the 
diagonals  do  not  depend  on  i  or  j. 

The  other  case  of  special  interest  is  the  first  bin  in  front  of  the  source.  This  gives  one  of  the 
diagonal  elements  of  the  matrix  W,  that  is,  the  element  Wjj; 

Wjj  =  1 .09  D  A"'-’  =  1 .09  D  '^  A  ( 15) 


where  Zj  =  iA.  The  formal  mathematical  problem  with  this  method  is  apparent  from  the  following. 
Consider  the  ratio  of  the  diagonal  element  to  an  off  diagonal  element  Wj^  where  j«i.  As  an 
example,  let  A  =  150  m  and  let  Zj  =  10  km.  Then 
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Wii/Wy  =  0.375  (A/Zj)’/3  =  9  X  lO"*. 


(16) 


The  peculiar  property  of  the  matrix  is  that  a  diagonal  element  is  smaller  than  any  other  element  in 
its  row.  The  diagonal  elements  can  be  smaller  by  several  orders  of  magnitude  This  lesults  from 
the  physical  path  weighting,  whereby  the  path  is  weighted  the  least  in  front  of  the  source,  zero  at 
the  source  and  increases  with  distance  away  from  the  source. 

The  numerical  problems  in  this  method  readily  exhibit  themselves.  Atmospheric  turbulence 
(and  hence  phase  distortions)  varies  by  four  orders  of  magnitude  from  a  C„^  value  of  about  10  '^ 
near  the  ground  to  10  '’  in  the  troposphere.  Optical  turbulence  can  vary  by  as  much  as  3  orders  of 
magnitude  in  the  troposphere,  depending  on  whether  there  is  a  strong  turbulent  layer  present.  The 
combination  of  the  range  of  turbulence  with  the  wide  range  in  the  matrix  W  requires  a  dynamic 
range  of  more  than  10’ m  the  receiver  detector  to  avoid  the  numerical  problems.  Alternatively,  the 
dynamic  range  requirements  can  be  reduced  by  limiting  the  resolution. 

Although  this  calculation  has  been  perfomied  for  a  specific  phase  effect,  similar  results  will 
likely  arise  for  other  phase  effects  as  well.  The  problems  arise  from  both  the  path  weighting  prob¬ 
lem  and  the  sensitivity  of  phase  effects  to  large  eddies.  These  two  effects  are  closely  related.  An 
alternative  approach  that  avoids  the  dominance  of  large  eddies  involves  those  effects  based  on 
scintillation.  For  a  point  source  at  altitude  Zj,  the  variance  of  log-amplitude  observed  by  a  receiver 
on  the  ground  is  given  by; 

a^2(Zi)  =  0.56  J^'  c„2(ti)  (q/Zi)^/^(Zj-q)-‘‘/''  dq.  ( 1 7) 

Note  that  the  weighting  is  zero  both  at  the  receiver  and  at  the  source  and  is  symmetric  about  the 
path  midpoint.  Exactly  as  in  the  phase  approach,  we  can  construct  a  set  of  triangular  equations 
corresponding  to  a  range-resolved  lidar  measurement.  The  resulting  matrix  W  gives  rise  to  the 
same  numerical  issues  as  the  phase  example. 

A  further  complication  arises  for  both  these  first-order  approaches  if  the  uplink  effects  are 
Included,  that  is,  the  propagation  from  the  laser  to  the  backscattering  region.  In  the  case  of  the 
phase  measurement  of  angle  of  arrival,  the  uplink  propagation  produces  beam  wander  at  the  back- 
scattering  altitude  Zj.  For  a  collimated  transmitter  of  diameter  a,  the  mean  square  wander  at  Zj  is 
given  for  propagation  from  z=0  to  Zj  by’ 

<pc^>  =  6.08  Zj^  a''^^  I  Cn^(h)dh.  (18) 

•'o 

The  illumination  of  the  backscattering  region  is  changing  as  a  result  of  the  beam  centroid  wander. 
The  image  of  the  backscattered  spot  wanders  in  the  image  plane  as  a  result  of  the  uplink  beam 
wander  and  the  downlink  image  dancing.  Clearly,  the  uplink  propagation  modifies  the  path 
weighting  function.  The  modification  is  not  a  straightforward  sum  of  the  uplink  and  downlink;  the 
up  and  down  paths  are  nearly  identical,  and  strongly  statistically  dependent.  Since  turbulence  can 
be  regarded  as  frozen  over  the  propagation  time,  the  entire  propagation  occurs  within  an 
isoplanatic  patch.  The  downlink  performs,  to  some  degree,  a  phase  conjugation  of  the  uplink  phase 
distortions,  thereby  reducing  the  dancing.  If  instead  of  backscattering,  a  retroreflector  is  used  at 


^Snkla,  R.J.  ( 1 988)  Vnifitd  Approach  to  Elerlromagnclic  Wave  Propagation  in  Turbulence  and  the  Evaluation  of  MuUiparameter  Integrals. 
Technical  Report  No.  807,  MIT/Lincoln  Laboratory. 
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altitude  Zj,  the  down  beam  would  be  a  phase  conjugated  version  of  the  up  beam  and  nearly  perfect 
conjugation  (except  for  diffraction)  would  be  observed  at  the  detector  with  almost  no  wander.  For 
backscattering,  the  down  beam  is  a  spherical  wave  (or  incoherent  superposition  of  spherical 
waves)  while  the  up  beam  is  a  collimated,  partially  coherent  beam  wave.  The  conjugation  is  less 
than  perfect.  This  is  a  difficult  problem  whose  solution  has  not  been  achieved.  It  is  clear  that  the 
path  weighting  function  has  been  modified  to  a  considerable  extent  by  the  inclusion  of  the  uplink 
propagation. 

If  this  overall  path  weighting  function  were  known,  it  is  still  likely  that  certain  difficulties  in 
a  phase  approach  would  remain.  The  problem  of  greatest  concern  is  the  requirement  of  very  large 
dynamic  range.  As  the  above  equations  show,  angle  of  arrival  fluctuations  are  directly  propor¬ 
tional  to  Cn^;  since  C„^  varies  by  more  than  10^,  so  do  phase  effects.  In  addition,  phase  effects  are 
dominated  by  large  turbulence,  which  occurs  near  the  aperture  and  in  the  boundary  layer. 
Measurements  from  backscattering  above  the  boundary  layer  will  include  these  large  distortions, 
and  the  contribution  from  above  the  boundary  layer  will  be  swamped  by  the  larger  contribution 
from  below. 

Closely  related  to  these  problems  are  the  noise  sources  inherent  in  a  phase  measurement  ap¬ 
proach.  One  source  is  laser  instabilities  that  contribute  a  jitter  to  the  propagation.  A  second  source 
is  mechanical  vibration  in  the  transmitter  or  receiver  that  also  contributes  a  jitter  and  results  in 
some  image  dancing.  In  the  presence  of  noise,  the  numerical  considerations  of  stability  become  of 
great  importance. 

The  value  of  scintillation-based  approaches  is  that  some  of  the  problems  inherent  in  the 
phase  approach  are  not  significant.  The  uplink  wander  of  the  beam  does  not  matter,  provided  the 
spot  remains  in  the  detector  field  of  view.  The  same  considerations  apply  to  other  sources  of  jitter, 
such  as  laser  instabilities  and  vibrations.  The  central  concept  in  log-amplitude  based  approaches  is 
that  they  are  insensitive  to  the  large  scale  sizes  and  the  important  scale  is  the  Fresnel  scale,  given 
by  (Xz)'^,  which  is  typically  of  the  order  of  centimeters.  Scintillation  along  an  optical  path  is  not 
dominated  by  the  largest  turbulence,  but  by  a  combination  of  the  value  of  C„^  and  the  turbulent 
layer’s  Fresnel  scale. 

A  key  issue  in  scintillation  is  that  the  conjugation  considerations  of  the  up  and  down  links 
are  a  negligible  consideration.  This  is  not  to  imply  that  scintillation  approaches  are  not  affected  by 
the  uplink  propagation.  The  uplink  propagation  results  in  intensity  fluctuations  of  the  light  inci¬ 
dent  on  the  backscatterers.  There  are  two  aspects  to  this  effect.  The  first  is  the  fluctuation  of  the 
total  incident  energy  on  the  backscatterers.  Scintillation  approaches  tend  to  calculate  the  mean 
square  average  of  the  log  intensity  or  log  amplitude  and  for  weak  scintillation,  these  can  be  shown 
to  be  virtually  the  same  as  the  normalized  intensity  variance.  The  approach  performs  a  normaliza¬ 
tion  for  the  overall  intensity  fluctuation.  The  second  aspect  pertains  to  scintillation  within  the  fi¬ 
nite  spot  size  of  the  beam  at  the  scatterers;  the  uplink  scintillation  will  result  in  intensity  hot  spots 
at  the  scatterers.  The  problem  is  that  of  scintillation  from  an  extended  source.  This  problem  is  a 
second  order  effect  and  is  discussed  in  the  following.  Closely  related  is  the  important  problem  of 
the  beam  broadening  due  to  the  uplink  distortions.  This  reduces  the  overall  scintillation  and 
changes  the  path  weighting.  Indeed,  the  use  of  different  size  sources  and  receivers  is  a  fundamen¬ 
tal  design  principle  in  the  horizontal  scintillometer  and  its  variants'".  Certain  effects  resulting  from 
beam  broadening  can  be  described.  The  path  weighting  function  for  scintillation  is  a  broad  peak. 


''’Och*,  O.R.  and  Wang,  T.  (1978)  Finite  aperture  acintillomeler  for  proniing  wind  and  Appl.  Opt.,  17:.^774-3778. 


8 


with  zeroes  at  the  transmitter  and  receiver.  The  location  of  the  center  of  the  peak  may  be  changed 
by  adjustment  of  the  receiver  and  transmitter  diameters,  but  the  peaks  will  remain  fairly  broad  and 
the  overall  scintillation  will  be  reduced  from  the  point  source  and  receiver  case.  Furthermore,  a 
practical  lidar  system  places  severe  constraints  and  limitations  on  the  amount  of  adjustment  pos¬ 
sible  in  the  receiver  and  transmitter  diameters.  These  considerations  will  result  in  the  previously 
stated  numerical  issues:  a  very  large  dynamic  range  in  the  receiver  is  required.  Very  limited 
altitude  resolution  would  be  achievable. 


4.  THE  CROSSED  PATH  TECHNIQUE 

The  ideal  lidar  measurement  method  would  be  based  on  an  optical  property  that  has  a  very 
sharply  peaked  path  weighting  function  and  would  be  minimally  effected  by  laser  instabilities, 
system  vibrations  and  uplink  beam  wander.  It  would  also  measure  an  optical  effect  that  yields  a 
reasonably  strong  signal.  The  first  requirement  insures  numerical  robustness.  It  has  also  been  sug¬ 
gested  in  the  previous  section  that  a  phase-based  approach  would  not  meet  any  of  these  require¬ 
ments  and  even  if  the  phase  conjugation  theory  for  the  uplink/downlink  problem  were  completed, 
the  sensitivity  to  system  jitter  would  remain.  It  appears  that  a  scintiUation-based  approach  holds  a 
higher  degree  of  promise  of  success. 

An  approach  is  presented  that  uses  a  scintillation  method  and  produces  a  sharp  path  weight¬ 
ing  function.  This  approach  is  an  adaptation  to  lidar  of  the  binary-star  methods  for  passive  remote 
sensing^’".  The  method  was  also  suggested'^  as  a  means  of  remote  sensing  across  a  horizontal  dis¬ 
tance  and  a  variant  has  been  used  to  measure  turbulence  in  fluid  flow'^.  The  strength  of  the  meth¬ 
od  is  that  the  path  weighting  is  determined  by  geometric  factors.  In  Figure  1,  two  parallel  beams, 
separated  by  a  distance  So,  are  backscattered  at  an  altitude  z.  The  backscattered  spots  are  viewed 
by  two  receivers  of  diameter  D,  separated  by  S,  in  such  a  way  that  their  viewing  paths  cross.  If  the 
altitude  of  the  backscatterers  is  given  by  z,  geometric  considerations  show  that  the  two  observation 
paths  cross  at  an  altitude  h  given  by 

h  =  Sz/(S+So)  =  z/(l+So/S).  (19) 

This  simple  equation  demonstrates  the  relation  between  the  backscattering  altitude,  the  crossing 
altitude  and  the  receiver  and  scatterer  separations.  The  second  form  shows  that  it  is  the  ratio  of  the 
backscatter  and  receiver  separation  that  determine  the  crossing  altitude.  In  a  lidar  implementation 
based  on  this  approach,  there  is  considerable  flexibility  in  adjusting  these  parameters. 

Each  receiver  measures  the  fluctuating  intensity  from  one  of  the  sources  and  this  scintilla¬ 
tion  is  determined  by  the  optical  turbulence  along  the  entire  path.  However,  if  the  correlation  of 
the  intensities  at  the  two  receivers  is  measured,  it  is  determined  solely  by  the  turbulence  in  com- 


''Vemin,  J.  and  Pelon,J.  (1986)  Scidar/lidar  description  of  a  gravity  wave  and  associated  turbulence:  preliminary  results,  Appl.  Opt.,  25:2874- 
2877. 

'^Wang,  T.,  Clifford,  S.F.  and  Ochs,  O.R.  (1974)  Wind  and  Refractive-Turbulence  Sensing  Using  Crossed  Laser  Beams,  Appl.  Opt.,  13:2602- 
2608. 

**Fisher.  M.J..  and  Krause,  F.J.(  1967)  The  crossed-beam  correlation  technique,  J.  Fluid  Mech.,  28:70.8-7 1 7. 
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mon,  that  is,  the  turbulence  in  the  volume  of  ovetl^  of  the  viewing  gromelries.  It  is  clear  that  the 
path  weighting  function  for  this  correlation  is  sharply  peaked  at  the  center  of  the  volume  of  over¬ 
lap,  given  by  h,  and  is  essentially  zero  outside  of  this.  The  geometry  can  be  used  to  provide  insight 
into  the  parametric  dependencies  of  the  width  of  the  weighting  function  and  hence  the  altitude 
resolution  of  the  technique.  If  A|  and  are  respectively  the  lower  and  upper  limits  of  overltq)  as 
shown  in  Figure  1,  then 

A,=D(z-h)/(S+So-D), 

(20) 

Aj  =  D(z-h)/(S+So+D). 

Since  in  practice,  the  receiver  diameter  would  be  smaller  than  the  other  horizontal  dimensions,  the 
approximation  D  «  S,  Sq  shows 

A  »  A,  =  Aj  =  D(z-h)/(S+So).  (21 ) 


Source  1 


Source  2 


r-h 


A  2 

A, 


Receiver  1 


D 

Receiver  2 


Figure  1.  Geometry  of  crossed  path  technique.  The  viewing  paths  of  the  two  receivers  are  crossed 
as  shown.  The  two  o  Jtgoing  lidar  beams  are  not  crossed  and  are  shown  by  the  dotted  lines. 
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This  is  equivalent  to  assuming  that  the  spherical  spreading  of  the  light  is  negligible  over  the 
volume  of  overlq>,  and  thus  the  light  cones  from  the  sources  can  be  treated  as  cylinders  in  the 
crossing  region.  As  a  numerical  exanq)le,  consider  backscattering  from  a  10  km  altitude.  If  the 
source  and  receiver  separation  are  1  m,  the  crossing  altitude  is  S  km  and  a  10  cm  receiver  gives  a 
geometric  overlap  of  +/-  250  m.  The  vertical  resolution  is  500  m. 

It  should  be  emphasized  that  these  geometric  relations  are  primarily  heuristic  and  are  given 
to  demonstrate  the  parametric  dependencies  and  semi-quantitative  behavior  of  the  solution.  A 
more  thorough  and  formal  solution  based  on  the  Rytov  theory  is  presented  shortly.  The  geometric 
theory  provides  only  semi-quantitative  information  about  the  path  weighting  function  and  its 
width,  while  the  Rytov  theory  specifies  the  function  completely.  For  the  numerical  example  just 
given,  it  will  be  shown  that  the  Rytov  theory  yields  much  sharper  path  weighting,  and  hence  high¬ 
er  resolution,  than  the  geometric  theory. 


4.1  Theory  of  the  Crossed  Path  Technique 


This  section  demonstrates  the  key  relationship  between  the  intensity  covariance,  which  is  the 
measurable  quantity,  and  the  covariance  of  log  amplitude,  which  is  related  to  the  crossing  volume. 
Let  1 1  and  Ij  be  the  intensity  received  by  the  two  detectors  shown  in  Figure  1.  In  the  more  general 
case,  both  sources  are  in  the  field  of  view  of  the  two  receivers.  This  minor  complication  is  ac¬ 
counted  for  later.  In  the  present  derivation,  we  assume  that  each  receiver  sees  only  one  source  and 
that  the  observation  paths  cross  as  suggested  in  Figure  1.  The  covariance  between  the  two 
receivers  separated  by  S  is 


C,(S)  =  <(I,-<I,»(l2-<l2»> 

(22) 

=  <I,l2>  -  <I,><l2>, 

where  the  brackets  indicate  ensemble  averaging.  By  definition  of  log  amplitude,  x.  we  have  that 
I,(x,y)  =  Ioi(x,y)  (23) 


where  Io|(x,y)  is  the  intensity  at  (x,y)  in  the  absence  of  turbulence.  The  analogous  relation  exists 
for  l2(x',y').  Upon  substitution  into  the  covariance,  we  get 


C,(S)  =  Io,(x,y)Io2(x’,y'){<e^^^'^^’y^  X2(3^'y)l>  .  <e2Xi(*.y)><e2X2(5^’.y’)>}^  (24) 


where  S  is  the  distance  between  (x,y)  and  (x',y’).  There  is  considerable  theoretical  and  experimen¬ 
tal  evidence  that  in  the  weak  scintillation  case  X  a  normal  (Gaussian)  random  variable  from 
which  it  follows  that 


(25) 


where  is  the  variance  of  log  amplitude.  Since  both  Xi  and  X2  are  normal,  so  is  their  sum,  x+, 
and  we  have  that 
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<e2X^  =  e2<X*>  +  V., 

where 

<X+>  =  <X«>  +  <X2>. 


(26) 


(27) 


and 


=  +  +  (28) 

This  last  form  introduces  the  covariance  at  a  separation  S  of  the  log  amplitude  at  receivers  1  and  2, 
Cj^(S),  which  is  formally  defined  by 


Cx(S)  =  <(Xt  -  <Xi>)(X2  -  <X2>)> 


=  <XiX2>  -  <Xl><X2>- 


(29) 


Putting  these  results  together  yields 

=  j2<x,>+2<X;>+2o,V2<j4+4C,(S) 

Substituting  these  results  produces 

C,(S)  =  {.2<x,>+2<x,>+2<t,;+2o4+4C,(S)  _  ^2<x,>+2<Xj>+2<J,;+2<J,^ 

=  I«,Io2  .  i} 

=  <I|Xl2>{e''*'x<^*- 1}.  (31) 

Approximating  the  exponential  for  weak  scintillations  gives 

C,(S)»4<I,><l2>C^(S).  (32) 

This  is  a  key  equation  in  our  formulation  which  demonstrates  that  the  intensity  covariance  be¬ 
tween  the  two  receiver  outputs  is  reduced  to  the  problem  of  evaluating  the  log  amplitude 
covariance.  It  is  worth  emphasizing  that  this  derivation  requires  only  the  assumption  of  normality 
of  the  log  amplitude.  Fried'^  performs  the  same  derivation  but  also  includes  the  requirement  of 
conservation  of  energy.  Since  optical  turbulence  does  not  dissipate  the  optical  energy,  conserva¬ 
tion  of  energy  is  implemented  by  requiring  that  <1>  =  1^,  which  in  turn  implies  that  for  log  normal 
variables,  <x>  =  As  the  above  derivation  shows,  the  relation  between  the  covariances  of  in¬ 
tensity  and  log  amplitude  does  not  require  that  conservation  of  energy  be  utilized.  It  is  clear  that 
this  restriction  on  <%>  conflicts  with  the  Rytov  approximation  which  produces  <x>  =  0  and  the 
normality  of  %.  The  Rytov  theory  is  only  an  approximation  resulting  from  a  perturbation  expan- 


'Tried,  D.L.  ( 1967)  Aperture  Averaging  of  Scintilliriona,  /.  Opt.  Soc.  Am.,  57: 169- 1 73. 


12 


sion  while  the  requirement  of  conservation  of  energy  is  exact.  Mixing  different  orders  of  approxi¬ 
mation  in  theory  is  to  be  avoided  since  it  may  lead  to  contradictions.  Since  the  theory  of  the 
crossed  path  technique  is  to  be  developed  from  the  Rytov  theory  and  the  conservation  of  energy 
requirement  is  unnecessary,  only  those  assumptions  that  are  part  of  the  Rytov  approximation  are 
assumed. 

Before  turning  to  the  derivation  of  the  log  amplitude  covariance,  it  is  shown  that  the  theory 
is  readily  expanded  to  accommodate  the  complication  of  both  sources  in  the  field  of  view  of  both 
receivers.  Let  I,  and  I2  denote  the  total  intensity  in  receiver  1  and  2  respectively.  In  the  case  of  in¬ 
dependent  sources,  the  total  intensity  results  from  an  incoherent  superposition  of  the  individual 
fields: 


(33) 

I2  =  I21  hi’ 

where  1^  is  the  intensity  from  the  j'*'  source  at  the  i*''  receiver.  Substituting  and  performing  algebra 
yields  the  following  form  for  the  covariance  of  the  intensity  between  receivers  1  and  2; 

Ci(S)  =  <Iil2>  -  <I|><l2>  (34) 

~  1 1  1 1  [■^^I2l22^"^l2^'^22^] 

['*'^l2l2l^**'J|2^*^2l^]  [‘^1|l22^‘*^ll^'^l22^]- 

In  this  expression,  the  terms  have  been  grouped  suggestively.  The  first  term  is  the  covariance  of 
source  1  at  receivers  1  and  2.  This  is  the  same  as  the  covariance  of  source  one  for  two  points  in  the 
observation  plane  separated  by  S  (if  the  receivers  are  regarded  as  point  receivers).  Now,  the  in¬ 
tensity  covariance  of  a  turbulent  layer  goes  to  zero  for  separations  larger  than  the  Fresnel  scale'-'. 
The  receiver  separation,  S,  will  in  practice  be  of  the  order  of  a  meter,  while  the  Fresnel  scale  of 
any  turbulent  layer  will  be  of  the  order  of  centimeters.  Therefore,  it  is  to  be  expected  that  in  the 
lidar  formulation,  the  first  term  is  small.  Similarly,  the  second  term  is  the  covariance  of  source  2  at 
a  separation  S,  and  is  also  small.  The  third  term  is  the  crossed  path  term;  it  is  the  covariance  of 
source  1  at  receiver  2  with  source  2  at  receiver  1.  It  is  given  by  the  earlier  result: 

<I  12121  ^"<112^^121^  ~  ^^01^02  (35) 

The  last  term  is  the  covariance  of  source  1  at  receiver  1  with  source  2  at  receiver  2.  This  is  the  un¬ 
crossed  path,  and  represents  the  covariance  of  independent  sources.  Again,  provided  the  paths  are 
separated  by  more  than  a  Fresnel  scale,  this  term  is  also  very  small.  Thus,  in  the  lidar  implementa¬ 
tion,  the  additional  terms  resulting  from  both  sources  in  the  field  of  view  of  both  receivers  are  ex¬ 
pected  to  be  small.  The  smallness  of  these  terms  in  the  above  arguments  is  relative  to  their  respec¬ 
tive  maximums,  and  not  relative  to  the  magnitude  of  the  crossed  path  term.  Although  the  above 


’’ctiffoid,  S.F.  ( 1978)  The  Clnitical  Theory  of  Wave  Propagation  in  a  Turboknt  Medium,  Chapter  2,  in  Laser  Beam  Propagalian  in  the  Atmo¬ 
sphere,  J.W.  Strohbehn,  Ed.,  Springer-Verlag,  New  York. 
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arguments  indicate  that  in  a  lidar  implementation  it  does  not  matter  whether  bodi  sources  are  in 
the  field  of  view  of  the  receivers,  this  claim  will  be  addressed  in  a  more  quantitative  fashion  in  a 
subsequent  section.  This  robusmess  is  a  very  desirable  property.  Because  of  this  pn^rty,  all  the 
following  derivations  will  assume  that  only  one  source  is  seen  by  each  receiver.  In  contrast,  the 
binary  star  technique  uses  receiver  separation  of  the  same  order  or  smaller  than  the  Fresnel  scale, 
so  the  additional  covariance  terms  must  be  included. 


4.2  Rytov  Theory  for  the  Weighting  Function 

The  theory  of  the  crossed  path  technique  has  been  reduced  to  that  of  the  covariance  of  the 
log  amplitude  for  the  propagation.  The  Rytov  solution  to  the  wave  equation  is  written  in  the  gener¬ 
al  form'* 

\|r(r)  =  k2[27iEo(r)]-'  n,(r*)  E;,(r')  dr’,  (36) 

•'  Ir-r'l 


where  y  =  x  +  i<t>i  X  is  the  log  amplitude  and  <{>  is  the  phase  and  thus,  x  =  Re(t|r).  In  this  equation,  k 
is  the  optical  wavenumber  (2ttA),  r  is  the  vector  observation  point,  Eo(r')  is  the  source  amplitude 
at  r',  and  n|(r')  is  the  refractive  index  fluctuation  at  r'.  By  definition  of  the  refractive  index  fluc¬ 
tuations,  <n,(r')>  =  0,  from  which  it  follows  from  the  Rytov  rqjproximation  thar  <\|/>  =  <x>  =  <♦> 
=  0.  Thus,  the  covariance  of  the  log  amplitude  reduces  to 

Cx(S)  =  <XiX2>-  (37) 

If  the  distance  to  the  observation  point  is  greater  than  any  lateral  dimension  in  the  problem, 
the  Fresnel  approximation  can  be  applied  and  so  the  Rytov  solution  can  be  wrinen  in  general  as 

V(r)  =  k2[27tEo(r)]-'  f  e^(P  ~  P')^/[2(z-z )]  En(r’)  dr’.  (38) 

•'  (z-z’) 


where  p  is  a  vector  in  the  plane  transverse  to  the  direction  of  propagation,  which  is  taken  to  be  the 
z-axis.  The  source  is  specified  by  Eo(r),  the  source  amplitude  observed  at  vector  r.  For  a  point 
source,  located  at  vector  I’d. 


Eo(r)  =  Eo|^jl^o'.  (39) 

If  the  point  source  is  located  at  z  =  0  and  transverse  vector  p^,  the  Fresnel  approximation  can  also 
be  applied  to  the  source  to  give 


Eo(r) «  Eo  e^  giklp-pol^/2z 

Tjcz 


(40) 


V.I.  (1961)  IVinr  Propagation  in  aTurhulenl  Medium,  McOraw-Hill.  New  York. 
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Substitution  into  die  Rytov  solutitm  gives 

V(r)  =  (k2/2jt)  z  f  il^-p’P/2(z-z’)  ^iklp'-poP/2z' 

zXz-z’) 

After  some  tedious  algebraic  manipulation,  it  can  be  shown  that 
Ip-PolVz  +  Ip-p’l7(z-z’)  +  Ip’-PfllVz*  =  lp'-Yp-Kl-Y)Pol^/y(2-z’). 
where  we  have  introduced  y  =  z'/z.  Making  this  substitution  gives 
V(r)  =  (k%n>  f  .Mp'-Yp-(l-Y)P,l^/2y<»-r’)  , 

To  simplify  this  further,  the  spectral  expansion  for  the  refractive  index  fluctuations  is  introduced: 

n(p,z)  =  J  dv(K,z)  e‘*^***,  (44) 

where  K  is  the  two-dimensional  vector  spatial  wavenumber.  If  this  substitution  is  made,  the  fol¬ 
lowing  integral  over  p'  results: 


(41) 


(42) 


(43) 


I,  =  Je**^*P  ^iklp'-'yp-(l-Y)polV2Tf(z-z') 

(45) 

By  a  change  of  variable  to  pa  where  p*  =  p'-YP-(l-Y)Po>  becomes 

=  reiK*P-jiklp.|2/2Kr-i') 

(46) 

=  2ne“‘-t1P*<‘-7)Pol  J J,(Kp.)  p.dp. 

(47) 

(48) 

The  first  step  in  this  integration  converts  the  integral  over  vector  p  into  an  integral  over  scalar  p 
and  angle.  This  allows  the  introduction  of  the  zerodi  order  Bessel  function  via  its  integral  repre¬ 
sentation‘s.  The  next  step  is  performed  by  recognizing  that  the  integral  is  in  the  form  of  a  Fourier- 
Bessel  (or  Hankel)  transform  of  a  Gaussian".  Using  this  result,  we  have  the  general  Rytov  solu¬ 
tion  for  a  point  source  at  arbitrary  vector  position,  Pg,  in  the  plane  z  =  0  and  observed  at  position 
vector  r  =  (p,z): 


i|r(r)  =  ik 


Jdz'J  dv(K^')  ,-iKr-r')K^/2k 


(49) 


'^MMhewn,  J.  Hid  Walker,  R.L.  ( 1970)  Methods  of  Mathematical  Physics,  2"*  edn.,  W.A.  Benjamin,  Reading,  MA. 
"Bracewell,  R.N.  (1978)  The  Fourier  Transform  and  Its  Applications,  2"^  edn.,  McOraw-Hill,  N.Y. 
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where  the  real  part  provides  the  log  amplitude.  Note  that  a  point  source  at  the  origin,  Po  =  0,  gives 


the  expression  usually  stated  as  the  Rytov  spherical  wave  solution.  It  follows  from  this  that  the  log 
amplitude  is  given  by 

X(r)  =  (l/2)[v  +  Vl  (50) 

=  k  Jdz' J  dv(K,z')  sin[7(z-z')K2/2k].  (51 ) 

The  requirement  that  the  refractive  index  fluctuations  be  real  is  used  to  obtain  the  last  expression. 
Specifrcally,  the  realness  of  n  in  the  spectral  expansion  implies  that 

dv(K,z)  =  dv*(-K,z).  (52) 

The  Rytov  solution  for  each  receiver  can  be  written  as 

X,  =  x(r,)  =  k  Jdz’J  dv(K,z')  sin[7(z-z')KV2k],  (53) 

X2  =  Xd*:)  =  k  Jdz"J  dv(K’,z")  e‘*^‘tT^P2+(l-y)Poil  sin[y(z-z")K’2/2k].  (54) 

where  y=z"/z.  The  subscripts  are  chosen  to  emphasize  the  crossed  path;  source  2  is  measured  at 
receiver  1  and  source  1  at  receiver  2.  Taking  the  product  of  these  two  expressions  and  performing 
ensemble  averages  gives 

<XiX2>  =  dz'dz"/  <dv(K,z’)dv(K',z")>  (55) 


giK»[7p,+(l-Y)po2]  sin[Y(z-z')K2/2k]  e‘*^’f^P2+(l-y)Poil  sin[y(z-z")K’2/2k]. 

The  rest  of  this  derivation  follows  exactly  the  ordinary  spherical  wave  solution”  and  will  not  be 
explicitly  done  here.  The  result  is 

<XiX2>  =  2tik2  Jdz’  J  dK  0„(K,z')  sin2[7(z-z’)K2/2k],  (56) 

where  ^>„(K,z')  is  the  three-dimensional  turbulence  spectrum.  For  the  Kolmogorov  spectrum,  this 
is 

«I»„(K,z)  =  0.033  C„2(z)K-"/^.  (57) 

Finite  inner  and  outer  scale  effects  can  be  included  by  using  the  modified  von  Karman  spectrum 

0„(K,z)  =  0.033  C„\z)  [K^  +  exp(-K^/K„^),  (58) 

where  is  related  to  the  inner  scale,  l^,  by  =  5.92/l(,  and  where  K^,  is  related  to  the  outer  scale, 
Lfl,  by  Ko=27tA-o.  Note  that  we  have  implicitly  assumed  that  neither  the  inner  nor  outer  scales  vary 
with  altitude. 


''’lnhhnani,  A.  (1978)  Wave  Propagation  and  Scattering  in  Random  Media,  Vol.  2,  Actdemic  Pre**,  N.Y. 
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The  expression  for  the  log  amplitude  covariance  can  be  further  simpliHed  by  using  the  fact 
that  the  turbulence  spectrum  depends  only  on  the  magnitude  of  K.  Consider  the  integral: 

I2  =  J  dK  <I>„(K,z’)  e^*^*P*  sin2[7(z-2’)K2/2k],  (59) 

where  we  have  defmed  p*  =  Y(Pi-P2H(l‘Y)(Po2'Poil-  Th®  vector  K  is  a  2  dimensional  vector  in  the 
plane  transverse  to  the  propagation  direction.  Let  the  direction  of  the  vector  p.  deHne  the  angular 
coordinate  0.  We  have 

I2  =  Jde  J  K  dK  <I»„(K,z')  sin2[Y(z-z*)KV2k] 

=  27C  J  K  dK  <I>„(K,z')  Jo(Kp.)  sin2[Y(z-z')K2/2k] . 

The  log  amplitude  covariance  is 

<XiX2>  =  4Tc2k2  Jdz’  J  K  dK  <D„(K,z’)  Jo(K|Y(PrP2)+(l-Y)(Po2-Poi)|)  sin2[7(z-z’)K2/2k]  (62) 

This  expression  can  be  simplified  through  the  specification  of  the  vector  positions  of  the 
receivers  and  backscattering  sources  for  the  crossed  path  geometiy.  The  receivers  and  sources  are 
assumed  to  be  symmetrically  located  about  the  origin  of  the  transverse  plane.  Furthermore,  it  is  as¬ 
sumed  that  the  two  sources  and  the  two  receivers  are  coplanar.  Without  loss  of  generality,  we  can 
define  the  sources  and  receivers  to  lie  on  the  x  axis.  With  reference  to  Figure  1,  we  have: 


Pi  “*P2 

Poi  ”  ■P02 

P2  =S/2 

Po!  =  S„n.  (63) 

If  we  specify  the  turbulence  spectrum  as 

<I»„(K,z)  =  C„2(z)<I)o(K),  (64) 

we  have 

<XiX2>  =  47c2k2  Jdz'  C„2(z')  J  K  dK  <Do(K)  Jo(KIyS-(1-y)SoI)  sin2[7(l-Y)zK2/2k]  (65) 

=  Jdz’C„V)W(z',z),  (66) 

where  the  substitution  (z-z')  =  (1-Y)z  has  been  made.  W  is  the  path  weighting  function  defined  by 
W(z',z)  =  47iV  j  K  dK  <I>o(K)  Jo(KIyS-(1-y)SoI)  sin2[7(  l-Y)zK2/2k].  (67) 


The  behavior  of  the  weighting  function  can  be  inferred  from  the  characteristics  of  the  in¬ 
tegrand.  In  Figure  2  the  log  amplitude  spectrum,  F,  is  plotted  where 


(60) 

(61) 
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F(K)  =  K  4>o(K)  sin2[*y(l-Y)zK2/2k] 


(68) 


and  4>o(K)  is  given  by  the  modified  von  Kannan  spectnim.  As  is  evident  from  the  figure,  the  spec¬ 
trum  is  insensitive  to  large  scale  sizes,  a  desirable  feature  of  any  scintillation-based  effect.  The  sin^ 
term  serves  to  attenuate  the  contributions  from  the  large  scale  sizes.  Also  noteworthy  is  the  broad 
peak  in  the  inertial  range,  with  a  maximum  approximately  given  by  the  first  maximum  of  the  sin^, 
that  is, 

Kf  =  [jik/'y(l-Y)z]  =  2n  [2X7(1-y)z]  ‘'2.  (69) 

This  wavenumber  corresponds  to  the  Fresnel  scale,  1^,  for  the  spherical  wave  case.  Using  a 
wavelength  of  0.5  )im,  this  scale  yields  a  numerical  value  of  5  cm  for  z  =  10  km,  z'  =  5  km.  Above 
this  scale,  the  log  amplitude  spectrum  is  rapidly  attenuated  by  the  turbulence  spectrum  and  the 
largest  contribution  to  the  spectmm  comes  from  Fresnel  scale  turbulent  eddies.  A  smaller  contri¬ 
bution  comes  from  larger  scales  in  the  inertial  range,  but  there  is  negligiMe  contribution  from 
smaller  scales.  Above  this  scale  the  spectmm  has  very  rapid  oscillations  due  to  the  argument  of 
the  sine. 

The  crossed  path  effects  enter  via  the  Jq  term.  Since  the  log  amplitude  spectmm  is  a  positive 
quantity,  it  follows  that  the  integral  is  a  maximum  when  Jq  is  a  maximum,  that  is,  when  the  argu¬ 
ment  of  Jo  is  zero.  This  condition  is  satisfied  when 


Figure  2.  Log  amplitude  spectrum  for  the  crossed  path  geometry.  The  spectrum  is  calculated  using 
an  outer  scale  Lq  of  10  m,  an  iimer  scale  lo  of  1  cm  and  A.  =  0.5  pm.  The  backscattering  altitude  is 
10  km  and  the  point  sources  and  point  receivers  are  each  separated  by  1  m.  The  crossing  altitude 
that  results  is  5  km  and  the  corresponding  value  of  the  Fresnel  scale  If  is  5  cm. 
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TfS  =  (1“Y)So, 


(70) 


or,  solving  for  z',  the  maximum  is  achieved  when  the  path  position,  denoted  by  h,  is  given  by 

h  =  Sz/(S+So)  =  z/(l+So/S).  (71) 

This  expression  is  identical  to  the  geometrical  crossing  altitude,  as  shown  in  Figure  1.  The  weight¬ 
ing  function  is  a  maximum  at  the  crossing  altitude;  the  geometrical  theory  and  the  more  exact 
treatment  agree.  Note  that  the  location  of  the  peak  depends  only  on  the  ratio  of  the  receiver  separa¬ 
tion  and  backscatterer  source  separation.  The  crossing  occurs  at  the  path  midpoint  if  these  are 
equally  separated. 

The  full  weighting  function  is  plotted  in  Figure  3  for  backscattering  from  10  km  and  equal 
source  and  receiver  separation  given  by  1  meter.  It  should  be  emphasized  that  this  is  the  point 
receiver  and  point  backscatterer  case.  The  weighting  function  is  very  sharply  peaked  at  5  km  and 
appears  to  be  negligible  at  less  than  4.5  km  and  more  than  5.5  km.  In  Figure  4,  the  vicinity  of  the 
peak  is  shown  on  an  expanded  scale  and  the  absolute  value  of  the  weighting  function  is  plotted. 
From  this  figure,  it  is  clear  that  the  weighting  function  has  a  1/e  width  of  the  order  of  175  meters. 
The  weighting  function  also  exhibits  a  very  rapid  decay;  beyond  about  300  m  from  the  peak,  the 
magnitude  is  down  by  nearly  a  factor  of  100  or  more,  and  beyond  800  m,  it  is  down  by  about  1000 
or  more.  This  width  does  not  depend  solely  on  the  ratio  of  the  receiver  to  source  separations,  but 
also  on  their  numerical  values  as  well.  It  would  be  satisfying  if  the  geometric  width  of  the  weight¬ 
ing  function  peak  also  agreed  quantitatively  with  the  exact  theory.  The  weighting  function  width 
depends  in  a  complex  way  on  the  first  zeroes  of  the  Bessel  function  and  no  simple  analytical  ex¬ 
pression  for  the  width  can  be  obtained.  A  more  thorough  discussion  of  the  resolution  of  the 
crossed  path  method  is  made  after  the  theory  is  extended  by  the  inclusion  of  finite  source  and 
receiver  effects. 

The  approach  of  the  crossed  path  technique  can  now  be  summarized.  Utilizing  our  earlier 
results,  it  follows  that 

[4<I|><l2>]''  QCS)  =  <XiX2^ 

=  Jdz'C„V)W(z’,z)  (72) 

The  intensities  are  readily  measured  by  photomultiplier  tubes  on  the  receivers.  The  average  in¬ 
tensities  and  their  covariance  are  statistically  estimated  from  the  time  series  of  each  receiver’s  in¬ 
tensity.  The  measurements  of  the  average  intensities  and  their  covariance  is  equivalent  to  the 
measurement  of  a  weighted  average  of  C„^  along  the  path.  These  characteristics  are  not  unique  to 
the  crossed  path  technique;  the  primary  feature  of  the  crossed  path  technique  is  the  sharpness  of  its 
path  weighting  function.  As  our  numerical  example  shows,  backscattering  from  point  sources  at 
10  km  with  1  m  separation  with  point  detectors  at  1  m  separation  gives  a  measurement  of  the 
weighted  average  of  in  a  175  m  layer  centered  at  5  km. 
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Figure  3.  The  crossed  path  weighting  function  versus  altitude  for  point  backscatterers  and  point 
receivers.  The  same  numerical  values  as  in  Figure  2  are  used. 


Figure  4.  The  absolute  value  of  the  weighting  function  is  shown  in  the  vicinity  of  the  crossing 
altitude  of  S  km.  The  same  numerical  values  as  in  Figure  2  are  used. 
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4.3  Finite  Receiver  and  Source  Effects 


It  is  well  known  that  scintillation  effects  are  reduced  by  receiver  aperture  averaging.  In¬ 
tensity  fluctuations  are  also  less  from  an  extended  source  than  from  a  point  source.  The  extension 
of  the  weighting  function  to  include  finite  source  and  receiver  effects  is  essential  to  the  design  and 
theory  of  a  lidar  C„^  measurement.  First,  the  extension  is  made  for  the  finite  receiver  effects.  It  is 
assumed  that  both  receivers  are  identical,  circular  apertures  of  diameter  D,  or  radius  R=D/2.  The 
measured  quantity  at  each  receiver  is  the  integral  of  the  intensity  over  the  aperture 

Iji  =  J  I|(x.y)  dxdy,  (73) 

where  1^,  denotes  the  total  intensity  and  an  analogous  expression  holds  for  1^2-  This  is  related  to 
the  log  amplitude  by 

Iti  =  J  dx  dy  (74) 

=  I„,  J  dx  dy, 

where  we  have  assumed  that  the  intensity  in  the  absence  of  turbulence,  I^,,,  is  uniform  across  the 
aperture.  Proceeding  as  before,  we  have  that 

<Iti>  =  loi  J  dx  dy 

=  lo,  7iR2  <e^5Ci>^  (75) 

where  the  local  homogeneity  of  x  has  been  used,  that  is,  <X(x,y)>  does  not  depend  on  (x,y)  within 
the  aperture  and  can  be  pulled  out  from  the  integral.  The  intensity  covariance  satisfies 


Cj  (S)  —  '<It|It2^  ■ 

It 


(76) 


=  Ifl,Io2  J  /  dx  dy  dx'dy'  X2(x’y)]>  .  <e2[Xi(x,y)><gX2(x',y’)>J 

Using  the  properties  of  the  normal  distribution,  it  follows  that 


Cj  (S)  =  [<Ij|><lT2>/(7tR^)^]  JJ  dx  dy  dx’dy'  ^ 

=  [<Iti><It2>/(^R^)^]  JJ  dx  dy  dx'dy’  4  Cj^(x-x',y-y').  (77) 

Since  =  <XiX2^’  **  **  clear  that  the  critical  quantity  consists  of  the  integral  of  <XiX2>  die 
apertures.  This  is  denoted  by  <XiX2>t  ^d  is  given  by 

<XiX2>T  =  JJ^^  dx’dy'  <Xi(x.y)X2(x',y')>-  (78) 


Note  that  it  is  normalized  by  the  square  of  the  receiver  area.  Interchanging  the  order  of  ensemble 
averaging  and  integration  gives 
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<XiX2>T  =  <[  J  dx  dy  Xi(x.y)]  [ J  dx’dy*  X2(*'.y’)]  >•  (79) 

The  problem  is  reduced  to  evaluating  the  integrated  log  amplitude  at  the  q>emues.  Call  these 
quantities  Xit  X2T  (o**  receivers.  These  are  given  by  the  forms 

XiT  =  ' J  dx  dy  Xi(x.y).  (80) 

where  the  integral  is  over  a  circle  of  diameter  D  centered  at  the  transverse  vector  p,.  Substituting 
for  Xi  from  our  earlier  results  produces 

X,T(r,)  =  k/(7cR2)  Jdp,  Jdz'J  dv(K,z')  e‘*^*t7(Pi+Pr)+(l-Y)Po2J  sin[Y(l-Y)zKV2k].  (81) 
The  integral  over  p,  inside  the  aperture  can  be  evaluated.  Consider  the  integral 

13=  Jdp,e'*^*7Pr.  (82) 


Converting  from  a  vector  integral  to  an  integral  over  angle  and  the  magnitude  of  p  produces 
I3  =  27C  J  p,  dp,  Jo(Kp,) 

=  27cR2 


(83) 


where  the  integral  has  been  evaluated  using  a  property  of  Bessel  functions.  The  effect  of  the 
receiver  aperture  averaging  is  to  include  a  further  term  in  the  integrand  that  behaves  as 
J,(K'yR)/K'^.  Using  this  result  in  our  expression  for  Xu- 

X,T  =  2k  Jdz’J dv(K,z')  e*'^*lVPl+<'-'.  ’Pozl  siii[Kl-y)zK2/2k].  (84) 

Note  that  the  scaling  by  the  receiver  area  nR‘  has  been  cancelled  out  by  our  normalization.  In  the 
following  we  will  drop  the  subscript  T  on  <XiX2^-  T^ie  derivation  of  <XiX2^  proceeds  exactly  as 
before  with  the  inclusion  of  this  additional  term  in  the  integrand  and  the  result  is 

<XiX2>  =  Jdz’  C„2(z')  J  K  dK  %{K)  Jo(KIyS-(1-Y)SoI)  sin2[Y(l-Y)z2/2k] 

(85) 

This  in  turn  defines  a  new  path  weighting  function  by 

W(z’,z)  =  IbnV  J  K  dK  %iK)  Jo(KIyS-(1-y)S„I)  sin^  [70-Y)^^/2k]  (86) 


The  J,(x)/x  term  in  the  integrand  acts  as  a  low  pass  filter.  Spatial  wavenumbers  satisfying  K  » 
1/(YR)  are  attenuated.  The  receiver  s^rture  averaging  reduces  the  scintillation  and  the  covariance 
by  fdtering  or  averaging  out  the  contributions  from  scale  sizes  smaller  than  the  aperture.  A  more 
precise  criterion  can  be  stated  using  the  1/e  point  of  [J,(x)/x]^.  This  value,  1.9,  is  used  to  define  the 
passband  of  the  filter.  ecifically,  wavenumbers  less  than  1.9/')R  will  not  be  attenuated  or 
reduced  by  aperture  averaging. 

The  receiver  aperture  averaging  is  a  critical  system  parameter.  A  larger  aperture  is  desirable, 
as  it  increases  the  total  signal  received.  It  cannot  be  so  large  as  to  significantly  reduce  the  scintilla- 
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tion.  These  two  criteria  can  be  satisfied  if  the  Fresnel  scale  eddies  are  within  the  passband.  Specif¬ 
ically,  we  require 

27c  [2X'y(l-Y)z]-‘^  ^  (1.9)/yR  (87) 

be  satisfied  for  the  altitude,  z',  where  the  weighting  function  is  peaked.  Additional  attenuation  at 
altitudes  not  in  the  weighting  function  peak  docs  not  matter.  Alternatively,  this  equation  can  be  in¬ 
verted  to  impose  the  requirement  that  the  aperture  radius,  R,  must  satisfy 

R  ^  (1 .9/27tY)  [2A,y(1-y)z]  •'2.  (88) 

In  practice,  the  backscattering  will  not  be  from  a  point  but  from  an  area  with  non-zero  cross 
sectional  area  determined  by  the  lidar  beam  width  at  the  scattering  altitude.  The  beam  may  wander 
and  be  broadened  due  to  the  turbulence  on  the  uplink.  These  phase  effects  will  be  discussed  sub¬ 
sequently.  Here,  the  scintillation  effects  of  the  uplink  are  considered.  The  lidar  beam  at  the  back- 
scatter  altitude  is  assumed  to  have  a  diameter  d  (or  radius  a  =d/2).  The  scintillation  on  the  lidar 
beam  results  in  a  fine  scale  speckle  pattern  of  intensity  within  a  circle  of  radius  a.  The  intensity  is 
not  uniform  within  the  circle  of  radius  a,  but  rather  there  are  intensity  hot  spots  of  the  size  given 
by  the  Fresnel  scale.  If  the  averaging  time  is  longer  than  the  wind  clearing  time  for  the  largest 
Fresnel  scale  hot  spots,  the  speckle  stmcture  of  the  light  incident  on  the  backscattering  region  can 
be  ignored  and  the  backscattering  region  can  be  assumed  to  be  uniformly  illuminated.  This  criteria 
can  be  more  precisely  stated.  The  lidar  beam  is  assumed  to  be  collimated  and  locally  fqjproxi- 
mated  by  a  plane  wavefront.  If  V  is  the  wind  speed  transverse  to  the  beam  at  z,  the  scattering 
altitude,  the  averaging  time  T  must  satisfy 

T>(kzym,  (89) 


where  the  plane  wave  Fresnel  scale  has  been  used.  For  a  0.5  pm  beam  at  10  km,  this  requires  that 
the  time  be  longer  than  0.07  sec  for  a  1  m/s  wind  and  7  msec  for  a  10  m/s  wind.  These  criteria  are 
easily  satisfied. 

The  illumination  of  the  backscatterers  can  be  regarded  as  uniform.  In  addition,  the  scattering 
process  is  either  Rayleigh  scattering  from  air  molecules  or  Mie  scattering  from  aerosols.  In  either 
case,  the  scattering  process  is  incoherent.  The  net  result  is  that  the  backscattering  of  an  incident 
beam  of  radius  a  can  be  treated  as  a  uniform,  incoherent  source  of  radius  a.  The  details  of  this  cal¬ 
culation  are  straightforward  but  tedious.  Let  I(x,y,r)  be  the  intensity  observed  at  (x,y,z)  from  a 
point  source  located  at  vector  position  r  in  the  z  =  0  plane.  For  an  incoherent  source,  observed  at 
(x,y,z),  the  total  intensity  1^  is  given  by  the  incoherent  superposition  of  point  sources; 


J  I(x,y,r)  dr. 


(90) 


where  the  integration  is  over  the  source  distribution  volume.  The  same  steps  as  for  receiver  aper¬ 
ture  averaging  are  followed  to  derive  the  effect  of  the  finite  source.  The  net  result  is  the  inclusion 
of  an  additional  term  in  the  integrand  for  the  log  amplitude  covariance,  a  term  given  by 
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The  resulting  log  amplitude  covariance,  incorporating  the  effects  of  finite  receiver  and  source,  is 
thus 

<X,X2>  =  Jdl'  C„V)  /  K  dK  J„(KIyS-(1-y)S„I)  (92) 

This  defines  the  overall  path  weighting  function  by 

W(z'^)  =  MnV  Jk  dK  4>„(K)  I„(KlyS-(l-Y)S„l)  (93) 

The  effect  of  a  source  of  radius  a  is  the  same  as  that  of  a  receiver  of  the  same  radius,  but 
with  Y  replaced  by  l-y.  The  finite  extent  of  the  source  results  in  filtering  of  the  spatial  scales  con¬ 
tributing  to  the  scintillation  spectmm  and  thereby  reduces  the  scintillation.  In  order  to  keep  this 
reduction  to  a  minimum,  we  again  consider  the  1/e  point  of  the  [J,(x)/x]^  term.  This  provides  the 
specific  criterion  that 

2n  [2Xy(  1  -Y)z] ' '^^  <  ( 1 .9 )/( 1  -Y)a  (94) 

be  satisfied  for  all  altitudes  of  interest.  The  requirement  for  the  backscatter  spot  size  is 

a  ^  1  .9/[27c(  1-y)]  [2>.y( I  -Y)z]  (95) 

The  following  examples  are  treated  to  provide  some  feel  for  the  numerical  constraints  on 
the  receiver  and  backscatter  spot  size.  The  laser  wavelength  is  assumed  to  be  0.5  pm,  the  back¬ 
scatter  spot  at  10  km,  and  the  crossing  altitude  is  5  km.  In  this  case  neither  the  lidar  spot  at  10  km 
nor  the  receiver  radius  can  exceed  approximately  3  cm  without  reduction  in  scintillation.  It  can  be 
seen  that  the  allowable  sizes  increase  as  the  square  root  of  the  backscatter  altitude.  For  a  2  km 
backscatterer  and  a  crossing  altitude  of  1  km,  the  limiting  size  is  about  1  cm.  These  are  tight 
restrictions.  In  Figure  5,  the  behavior  of  the  weighting  function  in  the  vicinity  of  the  peak  is  shown 
for  various  aperture  sizes  for  the  usual  numerical  example.  It  can  be  seen  that  the  effect  of  aperture 
averaging  is  to  cause  both  a  decrease  in  the  magnitude  and  an  increase  in  the  width  of  the  peak. 
The  signal  strength  and  resolution  are  degraded. 
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Figure  5.  The  absolute  value  of  the  weighting  function  is  shown  in  the  vicinity  of  the  crossing 
altitude  of  5  km  for  several  values  of  aperture  diameters.  The  same  numerical  values  as  in  Figure  2 
are  used,  except  that  the  receiver  and  backscatter  diameters  are  assumed  equal  and  given  by  the 
values  shown.  Note  that  the  excursions  to  small  values  shown  in  Figure  4  are  absent  in  the  present 
figure.  This  absence  is  a  numerical  artifact  of  the  larger  step  size  in  z  used  in  the  calculation. 


4.4  Inner  Scale  Effects 

In  the  theory,  one  scale  size  has  emerged  that  is  of  critical  importance:  the  Fresnel  scale. 
However,  the  presentation  has  oversimplified  one  aspect  of  the  problem;  namely,  we  have  as¬ 
sumed  that  the  inner  scale  of  turbulence  is  1  cm,  a  value  less  than  the  Fresnel  scale  in  the  calcula¬ 
tions  presented.  In  practice,  the  inner  scale  is  not  a  constant,  but  is  itself  a  function  of  both  altitude 
and  Uubulence  strength.  The  distinction  between  the  iruier  scale  for  velocity  and  temperature  fluc¬ 
tuations  must  be  made  in  discussing  inner  scale.  In  all  that  follows,  we  will  assume  the  inner  scale 
for  temperature  fluctuations.  The  inner  scale  for  temperature,  and  hence  refractive  index,  fluctuta- 
tions  is  defined  by^ 


lo  =  7.4(DVe)'^ 


(96) 


R.J.,  and  aifford.  S.F.  (1978)  Modiried  spectrum  of  atmospheric  temperature  fluctualiom  and  its  application  to  optical  propagation.  1. 
Opt.  Soc.  Am.,  ««:892-899. 
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where  D  is  the  diffusion  coefficient  and  E  is  the  eddy  dissipation  rate.  Values  for  1<,  have  been 
reported  as  millimeters^'  in  the  first  few  meters  above  the  surface,  centimeters  in  the  first  few  hun¬ 
dred  meters  and  in  the  troposphere^^  and  tens  of  centimeters  in  the  stratosphere^^. 

In  Figure  6,  the  value  of  the  weighting  function  is  shown  as  a  function  of  lo  for  different 
receiver  and  source  diameters.  The  value  of  the  weighting  function  is  plotted  at  the  crossing 
altitude  (for  the  S  =  S„  case)  where  the  weighting  function  is  a  maximum.  The  dependence  on  in¬ 
ner  scale  is  computed  via  the  modified  von  Kaiman  spectrum.  The  crossing  altitude  is  again  5  km 
and  yields  a  Fresnel  scale  of  approximately  5  cm.  As  is  evident  from  the  figure,  the  value  of  the 
weighting  function  does  not  depend  on  inner  scale  provided  is  smaller  than  the  Fresnel  scale. 
(The  effect  of  inner  scale  has  a  lesser  dependence  on  source  and  receiver  diameters.)  If  the  Fresnel 
scale  is  smaller  than  the  inner  scale,  the  refractive  index  spectram  dominates  the  integrand  of  the 
weighting  function  at  the  small  scale  sizes  of  the  inertial  range,  with  a  negligible  contribution  from 
the  sin^  term.  The  inner  scale  is  then  the  critical  parameter  instead  of  the  Fresnel  scale. 

In  Figure  7,  the  effect  of  inner  scale  on  the  weighting  function  in  the  vicinity  of  the  crossing 
altitude  is  shown.  Besides  the  suppression  of  the  magnitude  of  the  peak  value,  it  is  seen  that  there 
is  also  a  significant  broadening  of  the  weighting  function  when  1^  is  greater  than  the  Fresnel  scale. 
The  width  of  the  weighting  function  is  increased  and  hence  the  resolution  of  the  crossed  path  tech¬ 
nique  is  degraded. 

The  contribution  of  inner  scale  effects  to  intensity  variance  and  covariance  has  been  noted 
before^  for  propagation  through  homogeneous  turbulence  (that  is,  where  is  a  constant  along 
the  path).  Its  effect  in  the  crossed  path  technique  is  to  be  expected,  but  has  not  been  mentioned  in 
any  of  the  previous  works  on  the  variants  of  the  crossed  path  technique’’ Indeed,  it  ap¬ 
pears  that  the  effect  of  inner  scale  discussed  here  also  applies  to  the  techniques  based  on  binary 
star  data;  it  is  a  consequence  of  the  covariance  of  log  amplitude.  The  contribution  of  the  inner 
scale  in  a  crossed  path  lidar  system  must  be  eliminated  or  reduced  if  a  meaningful  measurement  of 
Cn’  is  to  be  made;  if  it  is  not,  the  measurement  depends  on  two  unknown  parameters.  If,  and 
The  situation  is  mitigated  somewhat  by  a  relation  between  E  and  C„^.  Specifically,  it  can  be  shown 
that 
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Figure  6.  The  maximum  value  of  the  weighting  function  is  shown  versus  the  inner  scale,  Iq.  The 
calculations  were  made  using  the  modified  von  Karman  spectrum.  The  source  and  receiver 
diameters  were  assumed  equal  and  given  by  the  values  shown.  The  source  was  located  at  10  km, 
giving  the  crossing  altitude  and  weighting  function  maximum  at  5  km.  The  other  values  are  identi¬ 
cal  to  those  used  in  Figure  2. 


Figure  7.  The  absolute  value  of  the  weighting  function  is  shown  in  the  vicinity  of  the  crossing 
altitude  of  5  km  for  several  values  of  inner  scale,  Iq.  The  same  numerical  values  as  in  Figure  2  are 
used,  except  that  the  receiver  and  backscatterer  diameters  are  assiuned  equal  and  given  as  S  cm. 
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C„2  =  ae2/3f(TpdT/^) 


(97) 


where  a  is  a  constant  of  order  unity  and  f  is  a  function  of  meterological  (that  is,  non-fluctuating) 
parameters.  Upon  substitution  of  it  follows  that 

C„2  =  0(7.4)*/^  D2  f(T,P,dT/dz)  (99) 

Thus,  C„^  is  inversely  related  to  !«;  large  C„^  means  small  l^.  We  have  also  seen  that  for  small  l^, 
the  Fresnel  scale  effects  dominate  and  the  weighting  function  depends  very  little  on  the  inner 
scale.  Strong  layers  of  turbulence  are  likely  to  have  small  l^,  (of  the  order  of  centimeters  or  less)  in 
the  boundary  layer  and  troposphere  and  the  weighting  function  can  be  unbambiguously  calculated 
for  these  cases.  For  weak  layers,  the  weighting  function  depends  on  the  unknown  l^,.  If  the  inner 
scale  is  assumed  to  be  less  than  the  Fresnel  scale  for  all  altitudes  of  interest,  the  effect  is  to  over¬ 
estimate  the  weighting  function  for  weak  layers,  and  thereby  underestimate  Cn^.  As  can  be  seen 
from  Figure  6,  the  averaging  effect  of  the  finite  sources  and  receivers  tends  to  reduce  the  over¬ 
estimation  of  the  peak  of  the  weighting  function. 

This  effect  of  inner  scale  cannot  be  eliminated  and  appears  to  impose  a  fundamental  limita¬ 
tion  on  the  crossed  path  method.  The  exact  contribution  of  inner  scale  is  unknown;  the  calcula¬ 
tions  performed  above  have  established  that  there  is  a  contribution,  but  these  calculation  have  been 
based  on  the  von  Karman  spectrum.  This  spectrum  is  ad  hoc,  without  experimental  basis,  and  was 
introduced  by  Tatarskii*'^  for  mathematical  convenience.  It  does  not  model  the  experimentally  ob¬ 
served  bump  in  the  temperature  spectrum^”.  Calculations  performed  by  Hill^*  for  a  point  source 
and  point  detector  show  for  propagation  through  homogeneous  turbulence  that  the  intensity  vari¬ 
ance  has  a  significant  contribution  from  the  inner  scale  region  that  is  manifested  for  l^  greater  than 
approximately  0.5(Xz)'^^.  Hill’s  analysis  may  be  interpreted  as  showing  that  the  modified  von  Kar¬ 
man  spectrum  is  a  reasonable  approximation  to  the  true  spectrum  for  l^,  greater  than  about  3(\z)'^. 
For  the  region  0.2(Xz)'^  to  2(Xz)'^.  the  agreement  is  poor  due  to  the  "bump"  and  the  modified  von 
Karman  spectrum  underestimates  the  spectmm.  A  further  interpretation  of  HUl’s  analysis  is  that  it 
provides  evidence  that  the  effect  of  inner  scale  may  be  more  prominent  than  the  effect  predicted 
by  using  the  modified  von  Karman  spectrum. 

The  criterion  for  the  inner  scale  effects  to  be  negligible  is  that  the  Fresnel  scale  be  larger 
than !(,,  or  specifically, 

1o<[2A.y(1-Y)z^]'^^  (99) 

where  z  is  the  altitude  of  the  backscattering.  The  contribution  of  the  inner  scale  of  the  spectrum  to 
the  crossed  path  covariance  may  be  minimized  by  maximizing  the  Fresnel  scale.  It  is  easy  to  see 
that  the  Fresnel  scale  is  maximized  for  y  =  1/2,  when  the  crossing  occurs  at  the  path  midpoint.  This 
is  equivalent  to  the  requirement  that  the  receiver  and  source  separations  be  equal.  This  provides  a 
strong  motivation  for  choosing  the  midpoint  as  the  crossing  point. 

Finally,  it  is  worth  noting  that,  in  the  limit  z  »  z',  the  above  requirement  then  becomes 

lo  <  (2Xz’)''2.  (100) 
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The  altitude  of  interest  is  now  the  crossing  altitude,  z',  rather  than  the  source  altitude.  This  is  the 
situation  with  the  binary  star  method.  Inner  scale  effects  should  manifest  themselves  in  the  binary 
star  data. 

It  can  be  argued  that  the  spatial  filtering  performed  by  the  finite  receiver  or  source  sizes  can 
be  used  to  filter  out  the  inner  scale  effect  on  the  covariance.  In  principle  it  is  possible  to  reduce 
these  effects  by  the  use  of  a  large  enough  aperture,  as  can  be  seen  from  Figure  6.  The  effect  will 
be  reduced  if  the  aperture  sizes  are  selected  such  that  the  spatial  wavenumbers  in  the  band  (27i/lo, 
Kf)  are  attenuated.  If  a  maximum  allowable  inner  scale  is  chosen  as  10  cm  and  equal  source  and 
receiver  separations  are  assumed,  the  filtering  of  these  effects  requires  receiver  and  source  sizes 
greater  than  about  12  cm  in  diameter.  A  practical  lidar  system  could  readily  meet  this  criteria,  but 
the  tradeoff  is  a  reduction  in  signal  strength  and  resolution.  In  addition,  the  filtering  provided  by 
J,(x)/x  is  gradual  and  does  not  have  a  sharp  transition;  therefore,  the  inner  scale  effect  is  not 
eliminated,  but  rather  reduced.  The  apertures  would  have  to  be  much  larger  to  get  a  significant 
reduction,  but  this  would  yield  a  significant  reduction  in  scintillation. 


4.5  Overall  Path  Weighting  Function 

The  complication  that  results  when  both  sources  are  in  both  receiver  fields  of  view  was  dis¬ 
cussed  in  Section  4. 1  in  terms  of  intensity.  The  weighting  function  described  so  far  has  applied  to 
the  crossed  path  term.  The  theory  can  be  extended  to  include  the  effect  of  the  other  three  terms  to 
yield  an  overall  path  weighting  function.  Equation  (34)  gives  the  intensity  covariance  between 
receivers  I  and  2; 

C[(S)  =  1 1 [^^I2^22^''^^I2^^^22^] 

+  [*~I|2l2l^*^^12^*^^2l^l  1 1^22^”^^ I l^*'^22^J *  (101) 

where  Ijj  is  the  j'*'  source  at  the  i'*’  receiver.  The  fonnal  expression  for  the  log  amplitude 
covariance.  Equation  (62),  is  in  fact  very  general  and  states  the  covariance  in  terms  of  arbitrary  2- 
dimensional  vectors  for  receiver  positions  p,  and  p,  and  for  source  positions  Poi  and  Poi-  The  first 
and  second  tenns  of  the  total  intensity  covariance  will  yield  a  path  weighting  function  found  by 
setting  Po,  =  p^j,  that  is,  both  sources  collocated: 

W,,2,(z',z)  =  W,2,22(z',z) 

=  647tV  J  K  dK  %(K)  Jo(KyS)  (102) 

where  the  subscripts  on  W  denote  the  particular  covariance  term  that  gives  rise  to  the  effect.  This 
path  weighting  represents  the  intensity  covariance  of  a  single  source  viewed  by  two  receivers 
separated  by  S.  The  third  term  is  simply  the  crossed  path  term.  In  our  current  notation,  this  term  is: 
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=  64)A2  J  K  dK  ♦oCK)  J„(KIiS-(1-7)SoI)  (103) 

I  »i>>^[T(l-r)iK"/2k] 

The  last  term  is  given  by  the  uncrossed  path,  which  we  get  by  setting  pg,  =  -pgj,  with  Po,  = 

So/2  (note  sign  reversal  from  crossed  path  case): 

W,  ,,22(z’,z)  =  64nV  J  K  dK  <Do(K)  Jo(KIyS+(1-y)SoI)  (104) 

[  ^  sin^[x  1  -y)zK2/2k] 

It  will  be  assumed  that  the  average  of  source  1  is  the  same  as  source  2  in  either  receiver  1  or  2. 
Specifically,  we  assume  that 

^I|  “  ^(21^  *  ^1 1 2^  *  ^l22^*  (10^1 

This  results  in  the  simple  addition  of  the  individual  path  weighting  functions  to  yield  the  overall 
path  weighting  function  given  by 

W(z',z)  =  64nV  J  K  dK  0„(K)  [2J„(KyS)  +  J„(KIyS-(1-y)S„I)  +  J„(KlYS+(l-y)S„l)  ] 

The  net  effect  of  both  sources  being  seen  by  both  receivers  is  the  inclusion  of  additional  Jg  terms 
in  the  integrand.  In  Figures  8  and  9,  the  overall  weighting  function  is  shown  for  comparison  with 
the  crossed  term  for  separations  of  10  and  1  cm,  respectively.  The  overall  weighting  function  for 
the  1  cm  case  is  significantly  different  than  the  cross  term  only.  For  the  10  cm  separation,  the 
weighting  function  at  the  peak  is  unaffected  by  the  additional  terms.  These  results  have  a  simple, 
physical  inteipretation  that  substantiates  the  qualitative  argument  presented  in  Section  4.1,  The 
relevant  Fresnel  scale  for  the  1  km  backscattering  with  midpoint  crossing  is  about  1 .6  cm  (for  X,  = 
0.5  pm)  corresponding  to  the  peak  in  the  log  amplitude  spectrum  (since  this  is  larger  than  the  inner 
scale  of  1  cm).  First,  consider  the  contribution  of  the  term  Jg(KYS).  If  Xg  is  the  first  zero  of  Jg,  there 
will  be  minimal  contribution  to  the  integral  over  K  for  values  satisfying 

K>VyS-  (107) 

Writing  this  in  terms  of  scale  sizes,  there  will  be  minimal  contribution  for  scale  sizes,  1,,  satisfying 

1,  <  27rYS/xg.  (108) 

Substitution  of  the  values  into  the  right  hand  side  for  the  peak  altitude  gives  1.3  cm  and  13  cm  for 
1  cm  and  10  cm  separations,  respectively.  The  spectral  peak  at  the  Fresnel  scale  is  within  the 
"passband"  of  the  Jg  term  for  the  1  cm  separation,  whUe  it  is  not  for  the  13  cm  case.  The  Jg  term 
that  arises  from  the  overall  path  scintillation  of  each  source  contributes  when  the  source  separation 
is  within  a  Fresnel  scale  and  does  not  contribute  significantly  at  larger  separation.  Note 
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Total  weighting  function 
Crossed  poth  term  only 
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Figure  8.  The  absolute  value  of  the  weighting  function  is  shown  in  the  vicinity  of  the  crossing 
altitude  of  SOO  m  showing  the  magnitude  of  the  total  function  and  the  crossed  term  only.  The  cal¬ 
culation  is  done  for  1  km  backscattering,  inner  scale  of  1  cm,  source  and  receiver  separations 
equal  and  given  by  10  cm,  and  source  and  receiver  apertures  of  1  cm. 


Figure  9.  The  absolute  value  of  the  weighting  function  is  shown  in  the  vicinity  of  the  crossing 
altitude  of  500  m  showing  the  magnitude  of  the  total  function  and  the  crossed  tenn  only.  The  cal¬ 
culation  is  done  using  the  same  numerical  values  as  Figure  8,  except  the  source  and  receiver 
separations  are  given  by  1  cm. 
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that  this  numerical  comparison  has  been  made  for  the  value  at  the  peak  corresponding  to  the  cross¬ 
ing  altitude.  Despite  the  negligible  contribution  at  the  peak  for  the  10  cm  case,  the  figure  shows 
that  for  other  altitudes,  the  contribution  from  the  uncrossed  terms  may  be  larger  that  that  from  the 
crossed  term,  but  these  are  less  important  since  the  weighting  function  is  down  by  more  than  an 
order  of  magnitude  from  its  peak  value.  The  contribution  of  the  JoCKyS)  term  is  not  symmetric 
about  the  midpoint,  but  applies  more  weight  to  the  receiver  end  of  the  path. 

Similar  arguments  apply  to  the  second  uncrossed  term,  Jo(KI7S+(1-y)SoI).  There  will  be  neg¬ 
ligible  contribution  from  this  term  from  scales  satisfying 

l2<27t[7S+(l-Y)So]/xo.  (109) 

Comparing  this  with  the  magnitude  of  the  term  from  the  single  source,  whole  path  scintillation,  it 
can  be  seen  that  1]  <  I2  for  any  values  of  S,  So,  and  Y-  From  this  inequality,  it  follows  that  this  sec¬ 
ond  uncrossed  term  has  a  narrower  passband  and  thus  a  smaller  contribution  than  the  single  source 
uncrossed  term.  In  our  specific  numerical  example,  the  values  of  Ij  are  2.6  cm  and  26  cm  respec¬ 
tively  for  1  and  10  cm  separations.  Both  these  values  are  larger  than  the  Fresnel  scale  and  thus  in 
both  cases  there  are  negligible  contributions  of  this  term  to  the  total  weighting  function. 

This  discussion  has  been  more  quantitative  than  the  earlier  arguments  about  the  importance 
of  the  uncrossed  path  terms  to  the  overall  weighting  function,  but  the  conclusions  have  been  the 
same.  The  earlier  discussion  argued  that  the  uncrossed  covariances  were  negligible  with  respect  to 
their  uncrossed  variances,  while  in  this  section  it  has  been  shown  that  they  may  also  be  negligible 
with  repect  to  the  crossed  covariance.  For  a  lidar  system,  the  separations  are  likely  to  be  tens  of 
centimeters  or  even  meters  '....c  ..tese  numbers  yield  cutoff  scale  sizes  for  the  uncrosssed  tenn  that 
are  larger  than  the  Fresn'*’  .‘■  •.Ue.  Consequently,  their  contribution  to  the  overall  weighting  function 
is  negligible.  This  is  not  true  for  the  binary  star  approach.  Indeed,  the  numerical  example  of  a  1  cm 
separation  at  1  km  is  equivalent  (or  10  cm  at  10  km)  to  an  angular  separation  of  the  sources  of  10 
prad.  This  is  the  separation  conesponding  to  the  double  star  Castor  which  has  been  used  in  the 
binary  star  method^'’. 


4.6  Uplink  Propagation  Effects 

A  critical  requirement  of  a  lidar-based  technique  that  has  not  been  addressed  is  the  effect  of 
the  uplink  propagation,  or  more  precisely,  the  effect  of  the  propagation  of  the  lidar  beam  from  the 
transmitter  to  the  backscattering  region.  Aspects  of  this  problem  have  been  mentioned  in  the  dis¬ 
cussion  of  the  effect  of  backscatter  spot  size.  It  was  argued  that  the  scintillation  effect  on  the 
uplink  would  not  be  a  problem  because  of  the  statistical  averaging  requirements.  With  respect  to 
beam  size,  a  very  strict  criterion  was  presented  in  terms  of  the  relevant  Fresnel  scale:  for  a  10  km 
backscatter,  the  (0.5  pm)  beam  should  be  less  than  6  cm  in  diameter,  while  for  1  km  backscatter,  it 
should  be  less  than  2  cm. 

This  criterion  needs  to  be  examined  in  terms  of  the  sj)ecific  effects  of  beam  broadening  due 
to  diffraction,  turbulence  and  beam  wander.  The  average  beam  width,  p^,  for  a  beam  of  diameter 
D,  can  be  expressed  in  terms  of  the  diffraction-related  width,  Pj,  and  the  turbulence-related  width, 
Pt.by 

<p^^>  =  (D/2)2ll-z/f]2  +  Pj2  +  <p,2>,  ( 1 10) 
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where  f  is  the  distance  to  the  focus  (for  a  collimated  beam,  f  =  «»)  and  where  the  diffraction  width 
is  given  by 

Pd  =  2z/(kD).  (Ill) 

The  turbulence  term  is  written  as^ 

<p  2>  =  4zV(kpo)2[l  -0.62(Po/D)'/^]«/5  (112) 

=  Pd'  (D/Po)"  [l  -0.62(Po/D)'/']«/^ 

where  Po  is  the  transverse  coherence  length-^  for  propagation  from  0  to  z  given  by 

Po=  |1.46k2  r C„2(t])(1  -Ti/z)5/3dnl-'/^  (113) 

•'o 

The  transverse  coherence  length  P(,  is  simply  related  to  the  atmospheric  coherence  length^'  r,,  by  r^ 
=  2.1  P(,.  This  expression  is  the  spherical  wave  fomi  of  the  transverse  coherence  length;  some  au¬ 
thors’  use  the  plane  wave  form  of  Po  for  a  collimated  beam  which  has  no  altitude  weighting  in  the 
integrand.  In  the  turbulence  literature,  this  beam  broadening  is  usually  discussed  in  terms  of  the 
short-term  beam  width  or  the  tilt-corrected  beam  width;  the  long-term  beam  width  includes  the 
turbulence  beam  wander.  Measurements  of  ro  for  the  plane  wave  case  yield  typical  values  at  0.5 
pm  of  from  5  to  10  cm,  with  higher  values  at  astronomical  sites.  These  values  translate  into  a 
range  for  p^  of  about  2.5  to  5  cm.  The  broadening  of  a  0.5  pm  beam  can  be  estimated  at  1  km  and 
at  10  km  using  the  value  of  2.5  cm  for  po-  This  requires  the  assumption  that  p^  is  the  same  for 
these  cases  as  for  the  ground  to  space  case.  This  is  reasonable  since  generally  the  largest  tur¬ 
bulence  occurs  near  the  ground  and  in  the  boundary  layer  and  dominates  Pp.  For  a  collimated 
beam  of  2  cm  diameter  at  the  transmitter,  the  turbulence  beam  broadening  will  be  much  smaller 
than  the  diffraction  effect:  the  2  cm  diameter  beam  will  be  of  the  order  of  2.6  cm  diameter  at  1  km 
and  17.4  cm  in  diameter  at  10  km.  For  a  focused  beam,  the  2  cm  transmitter  has  a  spot  size  of  1.9 
cm  if  focused  at  1  km  and  17.2  cm  if  focused  at  10  km.  A  5  cm  transmitter  produces  an  interesting 
example  where  both  turbulence  and  diffraction  effects  contribute:  at  10  km,  the  collimated  and 
focused  cases  yield  spot  sizes  of  1 1.6  and  10.6  cm,  respectively.  If  the  5  cm  transmitter  is  focused 
at  1  km,  the  spot  size  is  1.06  cm,  while  the  collimated  case  yields  5.1  cm.  The  spread  due  to  dif¬ 
fraction  can  be  reduced  by  increasing  the  aperture,  but  the  turbulence  spread  remains  and  is 
asymptotically  independent  of  the  diameter. 

These  calculations  demonstrate  that  beam  spreading  is  an  important  issue  in  the  crossed  path 
technique.  For  a  collimated  beam,  the  narrow  beam  widths  required  to  meet  the  specifications  for 
the  spot  size  result  in  diffraction  broadening  being  the  dominant  effect.  If  a  focused  beam  is  used 
instead  with  a  larger  transmitter  diameter,  turbulence  broadening  becomes  the  dominant  effect  in 


^Fanle,  R.L.  (1975)  Electromagnclic  Beam  Propagation  in  a  Turbulent  Media.  Proc.  IEEE.  63:  1669-1692. 

^Yura,  H.T.  (1971)  Atmospheric  Turbulence  Induced  Laser  Beam  Spread,  Opi.,  10:2771-2773. 

*'Fried,  D.L.  ( 1966)  Optical  Resolution  Through  a  Randomly  Inhomogeneous  Medium,  /.  Opt.  Soc.  Am.,  56;  1 372- 1 379. 
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determining  spK)!  size.  When  turbulence  dominates,  it  should  be  emphasized  that  the  beam  width  is 
a  fluctuating  quantity,  even  during  "homogeneous"  turbulence  conditions,  and  the  width  is 
specified  above  in  terms  of  its  ensemble  value.  Physically,  the  effect  of  beam  spread  is  to  reduce 
the  peak  value  of  the  weighting  function  and  broaden  its  width.  The  effect  of  the  fluctuating  lidar 
spot  size  is  that  now  the  weighting  function  itself  is  a  fluctuating  quantity.  The  relevant  quantity  is 
the  average  weighting  function.  This  in  turn  requires  that  the  transverse  coherence  length  be 
known,  or  equivalently,  the  profile. 

The  calculation  of  the  average  weighting  function  needs  to  be  done  numerically.  This  cal¬ 
culation  will  not  be  performed  here  and  is  left  for  follow-on  research.  Certain  aspects  of  the  aver¬ 
age  can  be  estimated,  however.  The  situation  is  not  as  significant  for  the  crossed  path  approach  as 
it  may  seem.  As  can  be  seen  from  Figure  5,  even  for  a  very  broadened  beam  of  20  cm,  the  weight¬ 
ing  function  is  still  very  narrow  and  far  superior  to  other  techniques.  Furthermore,  the  small  spots 
have  a  weighting  function  with  the  greatest  magnitude,  while  for  large  spots,  the  weighting  func¬ 
tion  is  the  smallest.  Whatever  the  distribution  of  spot  size,  the  average  weighting  function  will 
weight  the  less  probable  large  spots  the  least  and  spots  smaller  than  the  average  the  most.  This 
problem  needs  to  be  addressed.  Specifically,  the  distribution  of  the  beam  spread  is  required  to  cal¬ 
culate  numerically  the  average  weighting  function.  It  may  be  that  this  average  weighting  function 
does  not  significantly  differ  from  that  derived  by  using  the  RMS  spread  calculated  from  the  above 
expressions.  Another  possible  scenario  derives  from  the  apparently  circular  argument  that  the 
measurement  of  via  the  crossed  path  method  requires  the  knowledge  of  the  beam  spread  which 
in  turn  requires  the  knowledge  of  the  integrated  C„^.  The  crossed  path  measurement  may  need  to 
be  made  in  conjunction  with  transverse  coherence  length  measurements  of  a  stellar  source.  These 
auxiliary  measurements  would  then  be  used  to  calculate  the  spot  size,  and  hence  the  weighting 
function,  to  enable  the  crossed  path  technique  to  extract  altitude-resolved  C„^.  This  is  an  important 
aspect  to  the  problem  that  needs  to  be  examined  further. 

The  other  important  effect  of  the  lidar  uplink  propagation  arises  from  the  beam  wander.  The 
effect  of  beam  wander  is  that  the  backscatter  source  separation,  S^,,  becomes  a  random  variable. 
The  crossing  altitude  becomes  a  random  variable,  as  does  the  weighting  function.  Wander  results 
primarily  from  eddies  of  the  same  order  as  the  beam  size.  If  the  two  uplink  beam  are  separated  by 
a  distance  large  with  respect  to  the  beam  diameters,  the  wander  in  the  two  beams  can  be  consider¬ 
ed  as  statistically  uncorrelated,  since  different  eddies  will  be  responsible  for  the  wander  in  each 
upluik  beam.  The  wander  may  be  correlated  at  some  time  lag  due  to  frozen  turbulence  effects,  but 
the  correlation  at  the  same  time  is  treated  as  absent.  The  mean  square  wander,  of  either  a 

collimated  or  focused  beam  can  be  calculated  from 

<p>=  2.97  L^po’/^D'/^).  (114) 

It  follows  that  the  beam  wander  has  no  wavelength  dependence,  since  Po*''^  behaves  as  k'^.  This  ex¬ 
pression  shows  that  wander  can  be  reduced  by  increasing  the  transmitter  diameter.  Strictly,  this 
form  applies  to  the  regime  where  D  >  P(,.  For  the  case  where  D  <  p^,  there  is  negligible  wander 
and  the  primary  effect  of  turbulence  is  beam  broadening.  Wander  can  also  be  reduced  by  using  a 
small  enough  aperture,  that  is,  small  relative  to  the  turbulence  scales.  Fante^’  shows  results  from 
numerical  calculations  for  the  entire  range  of  D  and  P(,  that  exhibit  this  behavior. 

As  a  numerical  example,  consider  a  0.5  pm  beam  of  5  cm  diameter  at  the  transmitter,  with 
receiver  separation  of  1  m  and  the  nominal  backscatter  spot  separation  of  1  m.  For  propagation  to 
10  km,  the  nominal  crossing  altitude  is  5  km.  If  the  turbulence  is  specified  in  terms  of  2.5  cm  for 


Po,  the  RMS  wander  of  a  single  beam  is  calculated  to  be  5  cm  in  radius.  The  two  uplink  lidar 
beams  are  assumed  to  wander  independently  with  RMS  values  of  5  cm,  and  the  separation  be¬ 
tween  the  two  sf>ots  fluctuates  about  the  average  of  1  m  with  an  RMS  value  of  7  cm.  This  fluctua¬ 
tion  in  separation  yields  a  fluctuation  in  the  crossing  altitude  about  the  5  km  average  by  ±175  m. 
Clearly,  this  fluctuating  crossing  altitude  will  result  in  an  average  weighting  function  of  reduced 
peak  magnitude  and  greater  width  than  a  non-fluctuating  separation.  The  average  weighting  func¬ 
tion  due  to  wander  can  be  calculated  from  the  statistical  distribution  of  wander.  This  will  not  be 
done  here  and  represents  an  important  problem  for  follow-up  investigation.  As  with  beam 
broadening,  the  average  weighting  function  due  to  wander  is  neede^'  o  determine  and  gener¬ 
ally  requires  knowledge  of  the  wander,  which  in  turn  requires  po  oi  equivalently,  C„^  along  the 
path.  Such  data  can  be  provided  by  simultaneous  measurements  of  p^. 

Beam  broadening  and  wander  of  the  uplink  lidar  beam  are  important  parameters  that  aie  re¬ 
quired  to  determine  the  path  weighting  function  and  in  the  crossed  path  lidar.  The  impact  of 
these  parameters  must  be  taken  into  account  in  the  implementation  of  a  system  and  various 
tradeoffs  made.  For  example,  for  a  narrow  collimated  beam  of  the  order  of  a  few  centimeters 
diameter  at  the  transmitter,  wander  will  be  less  than  the  broadening.  In  effect,  the  beam  can  be 
considered  to  wander  within  a  radius  given  by  the  broadening  term.  Furthermore,  the  diffraction 
broadening  dominates  the  turbulence  broadening  in  such  cases.  For  such  narrow  transmitters,  the 
weighting  functions  can  be  calculated  using  our  earlier  expressions  using  the  diffraction-limited 
spot  size;  the  uplink  turbulence  effects  of  wander  and  broadening  can  be  ignored.  This  is  a  strong 
argument  in  favor  of  narrow  beams.  However,  these  diffraction-broadened  beams  wUl  exceed  the 
spot  size  requirements  discussed  earlier.  Such  narrow  beams  would  result  in  a  signal  and  resolu¬ 
tion  degradation.  This  situation  is  an  example  of  the  tradeoffs  required  in  a  practical  system. 


5.  CONCLUSION 

This  report  has  demonstrated  that  the  crossed  path  technique  possesses  unique  features  that 
can  be  exploited  for  the  lidar  remote  sensing  of  with  high  spatial  and  temporal  resolution.  The 
uniqueness  of  this  approach  stems  from  the  sharpness  of  the  path  weighting  function.  The  ap¬ 
proach  is  conceptually  a  geometric  one  and  can  be  regarded  as  a  form  of  intensity  interferometry. 
Since  the  approach  measures  the  normalized  intensity  covariance,  it  is  less  sensitive  to  fluctuations 
in  laser  output  and  laser  instabilities. 

There  are  numerous  parameters,  both  systemic  and  atmospheric,  to  take  into  account  in  a 
system  design.  There  are  many  tradeoffs  to  be  made.  One  imp>ortant  class  of  parameters  pertains  to 
the  strength  of  the  backscattering,  and  hence  the  signal  strength.  This  class  involves  the  choice  of 
wavelength  and  can  be  assessed  in  terms  of  a  simulation  using  such  existing  lidar  simulation  codes 
as  BACKSCAT^^.  The  strength  of  the  backscatter  signal  falls  off  as  density  in  the  case  where 
Rayleigh  scattering  dominates.  This  in  turn  determines  the  integration  time  and  the  temporal 
resolution  of  the  method.  In  order  to  compensate  for  this  decrease  in  signal  strength,  the  signal 
needs  to  be  strengthened,  perhaps  by  decreasing  the  receiver  or  source  separations.  This  would 
degrade  resolution.  Ideally,  the  resolution  of  the  weighting  function  would  exactly  match  the 


^^Ouiven*.  N.R.,  Rafune,  S.E.,  Hummel,  J.R.,  and  Cheifete,  M.O.  (1988)  BACKSCAT  Lidar  Backscalirr  Simulation  User's  Manual  for  Version 
1.0,  Oeophyaica  Laboratory,  Technical  Report  APOL-TR-88.033 1 ,  ADA2I9487, 
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receiver  gating.  With  repect  to  the  lidar  formalism  presented  in  Section  3,  this  situation  cor¬ 
responds  to  a  diagonal  weighting  function  matrix.  If  design  tradeoffs  require  the  use  of  broader 
weighting  functions  so  that  the  matrix  is  no  longer  diagonal,  the  crossed  path  approach  is  in  no 
way  invalidated.  The  path  weighting  function  is  still  very  sharp  and  localized  and  results  in  a 
sparse  matrix  (that  is,  one  which  has  only  a  few  nonzero  elements  adjacent  to  the  diagonal)  which 
is  readily  inverted. 

The  outlining  of  the  various  restrictions  on  the  system  parameters  seems  to  present  many 
conflicts.  It  is  doubtful  that  a  single  pair  of  receivers  could  perform  an  altitude  profile.  Rather,  a 
single  pair  of  receivers  could  be  used  to  probe  a  limited  altitude  range,  say  3  or  4  km  to  7  or  8  km 
with  a  vertical  resolution  of  150  meters.  Other  pairs  would  be  required  to  probe  other  altitude 
ranges;  this  is  not  regarded  as  a  serious  drawback,  since  the  receivers  consist  of  a  telescope  of  a 
few  centuneters  diameter  and  a  photomultiplier  tube  and  are  consequently  relatively  simple  and 
inexpensive.  The  achievement  of  this  goal  of  performing  some  prototype  experiments  to  measure 
via  the  crossed  path  methoa  over  a  limited  altitude  range  is  the  test  of  the  teclmique.  Tliis 
would  demonstrate  the  capability  to  remotely  sense  with  high  spatial  and  temporal  resolution, 
capabilities  that  currently  do  not  exist. 
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